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1. Introduction 



QCD, the theory of strong interactions, is an important component of the Standard Model 
of elementary particle interactions. It is asymptotically free, which allows us to compute 
cross sections of elementary particle interactions at high energies as a perturbative expan- 
sion in the running strong coupling a s (fiR). However, the running coupling a s (HR) remains 
rather large at the energies relevant at recent and future colliders. In addition, to leading 
order in the perturbative expansion, the coupling varies sizeably with the choice of the 
(unphysical) renormalisation scale /j,r. In hadron-initiated processes, the situation above 
is worsened by the dependence of the cross section on the (also unphysical) factorisation 
scale fj,p, which separates the long-distance from the short-distance part of the strong in- 
teraction. Thus, a leading-order evaluation of the cross section yields rather unreliable 
predictions for most processes in the theory. To improve this situation, in the past 25 years 
the radiative corrections at the next-to-leading order (NLO) accuracy have been computed. 
These efforts have culminated, when process-independent methods were presented for com- 
puting QCD cross sections to NLO accuracy, namely the slicing |l|, Q , subtraction ||, [|, |5| 
and dipole subtraction || methods. In some cases, though, the NLO corrections were found 
to be disturbingly large, and/or the dependence on fj,R (and eventually fj,p) was found to 
be still sizeable, thus casting doubts on the applicability of the perturbative predictions. 
When the NLO corrections are found to be of the same order as the leading-order predic- 
tion, the only way to assess the reliability of QCD perturbation theory is the computation 
of the next-to-next-to-leading order (NNLO) corrections. 

In recent years severe efforts have been made to compute the NNLO corrections to 
the parton distribution functions @ and important basic processes, such as vector boson 
production ^, 10, 11] and Higgs production || 12, 13, 14] in hadron collisions and jet 



production in electron-positron annihilation |15|, 16, 17]. These computations evaluate also 
the phase space integrals in d dimensions, thus, do not follow the process- independent 
methods used to compute the NLO corrections. Presently it is not clear whether those 



techniques IS, 19, 20, 21] can be directly applied for processes with more complex final 
states. 



The more traditional approach relies on defining approximate cross sections which 
match the singular behaviour of the QCD cross sections in all the relevant unresolved limits. 



Various attempts were made in this direction in Refs. [22, 23, |4|, 25, 26, 27, 



In general, the definition of the approximate cross sections must rely on the single and 
double unresolved limits of the QCD squared matrix elements. Although the infrared 
limits of QCD matrix elements have been extensively studied [31, 32, 33, 34], the formulae 
presented in the literature do not lend themselves directly for devising the approximate 
cross sections for two reasons. The first problem is that the various single and double soft 
and/or collinear limits overlap in a very complicated way and the infrared factorization 
formulae have to be written in such forms that these overlaps can be disentangled so 
that double subtraction is avoided. The second problem is that even if the factorization 
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formulae are written such that double subtraction does not happen, the expressions cannot 
straightforwardly be used as subtraction formulae, because the momenta of the partons 
in the factorized matrix elements are unambiguously defined only in the strict soft and 
collinear limits. In order to define the approximate cross sections one also has to factorize 
the phase space of the unresolved partons such that the singular factors can be integrated 
and the remaining expressions can be combined with the virtual correction leading to cross 
sections which are finite and integrable in four dimensions. 

In this paper we present a solution to the first problem. Due to momentum conservation 
constraints the expressions presented here are not unique. Nevertheless, these expressions 
may be useful in making the second step leading to a general subtraction scheme for 
computing QCD cross sections at the NNLO accuracy. We outline such a scheme for 
processes without coloured partons in the initial state, but do not explicitly define the 
approximate cross sections, which we leave for later work. 

The paper is organised as follows: after setting the notation in Sect. [2], we review in 
Sect. |3| the singly-unresolved limits and the subtraction terms relevant to NLO cross sec- 
tions. In Sect.|4| we review the doubly-unresolved limits and introduce the corresponding 
subtraction term. In Sect. |5] we introduce the subtraction term that regularizes the squared 
matrix element in all the unresolved regions of the phase space relevant to NNLO com- 
putations. That subtraction term is the keystone of this paper. In Sects. |6] and |7| we 
derive the iterated singly-unresolved limits and the singly-unresolved limits of the doubly- 
unresolved factorization formulae, respectively. Those singly-unresolved limits are used in 
the construction of the several contributions to the subtraction term presented in Sect. 
In Sect.|| we sketch the proof of the validity of the subtraction term of Sect.|| as a correct 
regulator of the divergences that occur in all the unresolved regions of the phase space rel- 
evant to NNLO computations. In Sect.0 we outline a possible general subtraction method 



for computing NNLO corrections to jet cross sections. In Sect. 1C we draw our conclusions. 



2. Notation 



We consider matrix elements of processes with m + 2 coloured particles (partons) in the 
final-state. Any number of additional non-coloured particles is allowed, too, but they will 
be suppressed in the notation. Resolved partons will be labelled by k, . . ., unresolved 
ones are r and s. 

The colour indices of the partons are denoted by c«, which range over 1, . . . , N% — 1 for 
gluons (or any other partons, such as gluinos, in the adjoint representation of the gauge 
group) and over 1, . . . , N c for quarks and antiquarks (or any other partons, such as squarks, 
in the fundamental representation). Spin indices are generically denoted by Sj. As in 
Ref. || , we formally introduce an orthogonal basis of unit vectors \c\, . . . , c m ) (g) \s\, . . . , s m ) 
in the space of colour and spin, in such a way that an amplitude of a process involving 
m external partons, Mm USl ^ {{pi}) with definite colour, spin and momenta {pi} can be 
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written as 

Mt' Sl} ({Pi}) = ((ci, • • • ,Cm\ ® {si, ■ ■ ■ , s m \) \M m ({pi})) . (2.1) 

Thus \M m ) is an abstract vector in colour and spin space, and its normalization is fixed 
such that the squared amplitude summed over colours and spins is 

\M m \ 2 = (M m \\M m ) . (2.2) 

The matrix element has the following formal loop expansion: 

\M) = \M {0) ) + \M {1) ) + |M (2) ) + . . . , (2.3) 

where |.M( )} denotes the tree-level contribution, (.M^ 1 )} is the one-loop contribution, 
\A4^) is the two- loop expression and the dots stand for higher-loop contributions. The 
amplitude \M) is assumed to be renormalized. In this paper we study the infrared be- 
haviour of the tree-level contribution, therefore, renormalization concerns us only to the 
extent of regularization. We use tree amplitudes obtained in conventional dimensional 
regularization. 

Colour interactions at the QCD vertices are represented by associating colour charges 
Tj with the emission of a gluon from each parton i. The colour charge T, = {T™} is a 
vector with respect to the colour indices n of the emitted gluon and an SU(N C ) matrix 
with respect to the colour indices of the parton i. More precisely, for a final-state parton i 
the action onto the colour space is defined by 

(c 1 ,...,Ci,...,c m \Tf'\b 1 ,...,bi,...,b m ) =S clbl ...T£ b ....6 Cmbm , (2.4) 

where is the colour-charge matrix in the representation of the final-state particle i, i.e. 
T cb = i/cnfe if i is a gluon or a gluino, T™^ = t™p if i is a (s)quark and = -t^ a if i is 
an anti(s)quark. Using this notation, we define the two-parton colour-correlated squared 
tree-amplitudes, |-M(^ fc )| 2 ) as 

l-^(°i)(fe})N(-^ (0) ({w})l^-^!-^ (0) ({w})) 

J i % k K 

and similarly the doubly two-parton colour-correlated squared tree-amplitudes |A4? fc N ^| 2 , 

\ M mw/= ( M(0) \{ T i ■ Tk,Tj ■ Ti} \M<®) , (2.6) 

where the anticommutator {Tj-T^, Tj-Ti) is non-trivial only if i = j or k = I, see Eq. (|2.§|). 

In our notation, each vector \A4) is a colour-singlet state, so colour conservation is 
simply 

(j2T^)\M) = 0, (2.7) 
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where the sum over j extends over all the external partons of the state vector \M), and 
the equation is valid order by order in the loop expansion of Eq. (|2.3| ) . 

The colour-charge algebra for ^2 n T™T£ = Ti-T k is 

T i -T k = T k -T i if i + k- Tl = d. (2.8) 

Here C% is the quadratic Casimir operator in the representation of particle i and we have 
C F = T R (iV c 2 - 1)/N C = (A c 2 - 1)/(2N C ) in the fundamental and C A = 2T R iV c = N c in the 
adjoint representation, i.e. we are using the customary normalization Tr = 1/2. 



3. Factorization in the singly-unresolved collinear and soft limits 

As explained in the introduction in order to devise the approximate matrix elements, we 
have to study the factorization properties of the relevant squared matrix elements when 
one, or two partons become soft or collinear to another parton. The factorization of the 
tree-level matrix elements when one parton becomes soft or collinear to another one are 



well known [35, 36]. In order to introduce a new notation, in this section we recall the 
relevant formulae for a tree-level matrix element |A4^:l 2 ) with m + 2 massless QCD partons 
in the final state and two massless and colourless particles in the initial state. 

3.1 The collinear limit 

We define the collinear limit of two final-state momenta pi and p r with the help of an 
auxiliary light-like vector nf r (nf r = 0) in the usual way, 

1,2 fl ju2 n 

Pi-ZiPir K r±_ 9 » Pr - ZrPir + K r± 7 o_. _. ' \ 6 - V ) 

where pf r is a light-like momentum that points towards the collinear direction and k r ±_ is 
the momentum component that is orthogonal to both pi r and rij r (pi r ■ k r ± = rii r ■ k r ± = 0). 
Momentum conservation requires that z% + z r = 1. The two-particle invariant masses 
(si r = 2pip r ) of the final-state partons are 

k 2 L 

Sir = • (^"^) 

The collinear limit is defined by the uniform rescaling 

krj_ — \k r j_ , (3.3) 

and taking the limit A — > 0, when the squared matrix element of an [m + 2)-parton process 
has the following asymptotic form, 



\M^ +2 (p hPr , . . .)| 2 ~ 87ra s ^ i- (M { Zi(Pir, ■ ■ .)\P$ r (zi,z r , krr,e)\M { 2{Mr, ■ ■ ■)) 



(3.4) 



-6- 



In Eq. the meaning of the ~ sign is that we have neglected subleading terms (in 
this case those that are less singular than 1/A 2 ). The (m + l)-parton matrix elements 
on the right-hand side of Eq. fl3.4| ) are obtained from the (m + 2)-parton matrix elements 
by removing partons % and r and replacing them with a single parton denoted by ir. 
The parton ir carries the quantum numbers of the pair i + r in the collinear limit: its 
momentum is pf r and its other quantum numbers (flavour, colour) are obtained according 
to the following rule: anything + gluon gives anything and quark + antiquark gives gluon. 
The kernels Pj°J are the ci-dimensional Altarelli-Parisi splitting functions, which depend on 
the momentum fractions of the decay products and on the relative transverse momentum of 
the pair. For the sake of simplicity, we label the momentum fraction belonging to a certain 
parton flavour with the corresponding label of the squared matrix element, Zh = z%. In the 
case of splitting into a pair, only one momentum fraction is independent, yet, we find it 
more convenient to keep the functional dependence on both Zi and z r . Depending on the 
fi flavours of the splitting products the explicit functional forms are 

IV I th I 



(At|pW p (^,z r ,fcj.;e)|i/) = 
(r\Pi%(z t ,z r ;e)\ S ): 



2C A 



9' 



£1 + £1 

Z r Zj, 



2(1 — e)ziZ r 



A- 2 



gT + 4ZiZ r 



rv i A/ i 



S rs Cp 



1 + zf 



ez r 



(3.5) 
(3.6) 
(3.7) 



where in the last equation we introduced our notation for the spin-averaged splitting func- 
tion, 

Pfifrfa'Zrl 6 ) = { P fifr{Zi,Zr,kx]E)) . (3.8) 

The gluon-gluon and quark-antiquark splittings are symmetric in the momentum fractions 
of the two decay products, while the quark-gluon splitting is not. Nevertheless, we do 
not distinguish the flavour kernels P qg and P gq . The ordering of the flavour indices and 
arguments of the Altarelli-Parisi kernels has no meaning in our notation, i.e., 



P fifr( Z i> Z r\ £) = Pfrfi( Z r,Zf,£) . 



(3.9) 



Thus, it is sufficient to record the kernel belonging to one ordering. We keep this convention 
throughout. 

In order to simplify further discussion, we introduce a symbolic operator Cj r that 
performs the action of taking the collinear limit of the squared matrix element, keeping 



the leading singular term. Thus we can write Eq. (3.4) as 



C ir \M { V +2 ( Pi , Pr ,...) 



87ra: s /i 



2o 



— (M 

Sir 



(0) 

m+1 ( 



(p ir ,...)\P$\M% + l (Pir, ■■■))• (3-10) 



,(0) 



fif: 



3.2 The soft limit 



The soft limit is defined by parametrizing the soft momentum as p r = Xq r and letting 
A — > at fixed q r . Neglecting terms that are less singular than 1/A 2 , the soft limit of the 
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squared matrix element can be written as 

\M^ +2 {p r , ...)\ 2 ^ -8na s ^ £ \s ik {r)\M { ^ +Hhk) {. . .)| 2 , (3.11) 

i,k 

where 

S tk (r) = (3.12) 



is the eikonal factor. In Eq. ( 3.11[ ) the (m + l)-parton matrix element on the right-hand 



side is obtained from the {m + 2)-parton matrix element on the left-hand side by simply 
removing the soft parton. 

Similarly to the Cj r operator, we introduce the symbolic operator S r that performs the 
action of taking the soft limit of the squared matrix element, keeping the leading singular 
terms. Thus we can write Eq. ( |3.11 ) as 



Srl-M£U(Pn • • -)! 2 = -8na s ^ £ \s ik {r)\M^ +l .^ k) {. . .)| 2 , (3.13) 

i,k 

if r is a gluon and S r \Ai^ +2 (p r , . . .)| 2 = if r is a quark. 
3.3 Matching the singly-unresolved limits 

If we want to regularize the squared matrix elements in all singly-unresolved regions of the 
phase space then we have to subtract all possible collinear and soft limits, i.e. subtract the 
sum 

E (E^ C - + S -| \M^ +2 ( Pi , Pr ,...)\ 2 , (3.14) 

where the 5 symmetry factor takes into account that in the summation each collinear 
configuration is taken into account twice. Subtracting Eq. ( 0.14 ) we perform a double 



subtraction in some regions of the phase space where the soft and collinear limits overlap. 
The collinear limit of Eq. ( 3.13j ) when the soft gluon r becomes simultaneously collinear to 
parton i is 

C ir S r \M^ +2 (p h p r , • • Ol 2 = -8vra s ^— £ J|A^> k) ( Pl , . . .)| 2 . (3.15) 

Sir ^ Z r 



The factor z%jz r is independent of k, therefore, using colour conservation (Eq. (2/7)) we 
can perform the summation and obtain, 

C tr S r \M^ +2 (p t ,p r , . . .)| 2 = STra^—^TjlM^iiPi, • • Ol 2 • (3-16) 

Sir %r 



Similarly, the soft limit of Eq. ( 3.10 ) when r is a gluon and z r — ► is 



S r C ir \M^ +2 {pi,p r , • • .)| 2 = 87ra sf i 2£ — T 2 -\M^ +1 ( Pu ...)\ 2 . (3.17) 

S{r #r 
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Eqs. (|3~16l) and (|317l ) differ by the factor z { = 1 — z r in the numerator of Eq. ( |3.16[) , in 



which z r becomes subleading if r is soft. Therefore, Eq. ( 3.1 6| ) can be used to account for 
the double subtraction: it cancels the soft subtraction in the collinear limit by construction, 



m+2l 



o, 



and the C? r — C^ r S r difference is subleading in the soft limit, 



(0) 1 2 
m+2l 



0. 



(3.18) 



(3.19) 



Accordingly, in order to remove the double subtraction from Eq. ( |3.14 ), we have to add 
terms like that in Eq. ( 3.16| ), That amounts to always take the collinear limit of the soft 
factorization formula rather than the reverse (like terms in Eq. ( |3,17| )). Thus the candidate 
for a subtraction term for regularizing the squared matrix element in all singly-unresolved 
limits is 



Ax\M 



m+2\ 



S r ^ Ci r S r 



\M<®. 2 (pi,Pr,...) 



iri ...)\P$ r \MZl(Pir,---)) 



(0) 



- S fr9 E 

k^i,r 



2_Zi_ 

Sir Zr 



l-^L+i ; (;,fc)(- 



(3.20) 



Note that the cancellation of the collinear terms in the soft limit actually requires the 
symmetry factor multiplying the collinear term, but not the collinear-soft one. This feature 
will be valid throughout the paper. 



The formula in Eq. ( 3.20 ) cannot be used as a true subtraction term. The reason is that 
the factorization formulae are valid in the strict collinear and/or soft limits. In order to 
define subtraction formulae over a finite part of the phase space, we either have to specify 
the momentum that becomes unresolved (general solutions are presented in Refs. ||, Q), or 
exact momentum conservation has to be implemented (as in Ref. ||). In the latter case in 
addition to the poles that are shown explicitly, these subtraction terms contain singularities 
in the (m + l)-parton matrix elements, 



C,sAi|A^ 2 | 2 = 



8,M\M 



where 



C ? ' S Ai 



m+2\ 



C JS )\M^\ 2 



2 + ( S S A 1 - S s 



\M {0) I 2 

\ Jvl m+2\ J 



(3.21) 
(3.22) 



^ Cj S ( S r C sr s r Cj, r s r + ^ ^ ^— Cir Cj r s r ^ I 



\M^ +2 1 2 , (3-23) 
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S s Ax 




1 2 - (3-24) 



The terms in Eqs. Q3.23Q and (ft.24Q are the iterated singly-unresolved limits that we shall 
define precisely in Sect. |6|. In a NLO calculation these terms give singularities in the 
doubly-unresolved parts of the phase space, where the measurement function, that defines 
the physical quantity (see Sect.|9|), becomes zero, and thus screens the divergencies. The 
same is not true in a NNLO calculation, which we shall tackle in Sect.|8|. 



4. Doubly-unresolved limits 



In this section we discuss all the various limits when partons r and s are simultaneously 
unresolved, which are needed for constructing an approximate cross section for regularizing 
the squared matrix element in the doubly-unresolved regions of the phase space. The 
infrared factorization in the various doubly- unresolved regions was discussed in Ref. p7| , 
but without aiming at combining the different factorization formulae such that double 
counting among the various expressions relevant for single and double soft and/or collinear 
limits is avoided. In order to set our notation and also exhibit these common terms, we 
rewrite the published formulae. 

At NNLO, the dependence of the squared matrix element |.M m +2| 2 on the momenta of 
the final-state partons after integration over the phase space leads to leading singularities 
in five different situations, 

1. the emission of a collinear parton-triplet; 

2. the emission of two pairs of collinear partons; 

3. the emission of a soft gluon and a pair of collinear partons; 

4. the emission of a soft quark-antiquark pair; 

5. the emission of two soft gluons. 

The first case is the triply-collinear one, when three final state partons become collinear. 
The next one corresponds to the doubly-collinear region, where two final-state partons 
become collinear to two other final-state partons, but not to one another. The third case 
occurs in the soft-collinear region, where the momentum of a gluon becomes soft in any 
fixed direction and, at the same time, two partons become collinear. The last two situations 
occur in the soft region, where either the momenta p r and p s of a final-state quark-antiquark 
pair tend to zero, with the directions of p r and p s as well as p r /Ps fixed, or the momenta 
of two final-state gluons tend to zero in any fixed direction. In the following we recall the 
most singular behaviour of the squared matrix element in those regions of phase space. 
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4.1 Emission of a collinear-parton triplet 



We consider three final-state partons i, r and s that can be produced in the splitting process 
irs — > i + r + s. The partons i, r and s have momenta pi, p r and p s . We introduce the 
following Sudakov parametrization of the parton momenta, 



2-j ^Pirs^ 



where p^, s is a light-like momentum that points towards the collinear direction, n^ rs is 
an auxiliary light-like vector (nf rs = 0) and kj± are the momentum components that are 
orthogonal to both pi rs and rii rs (pi rs ■ kj± = rii rs ■ kj± = 0). The variables Zj and kj± 
satisfy the constraints Zi + z r + z s = 1 and ki±_ + k r ± + k s x_ = 0. 

The two-particle invariant masses of the final-state partons are 



kj± h± \ 



2 



(Pj + Pi) = ~ z i z i H . j,l = hr,s. (4.2) 

The triply-collinear region is identified by performing the uniform rescaling 

k j± -> A%_l , j = i,r,s , (4.3) 

and studying the limit A — > 0. Again, we introduce the symbolic operator Cj rs that keeps 
the leading singular (0(1/A 4 )) terms of the squared matrix element in this limit. According 
to Ref. Wh i the squared matrix element |A4 m+ 2| 2 behaves as, 



Ci rs \M^ +2 {pi,Pr,Ps, ■ ■ OP 



(8ira s ^) 2 (MW( Ptrs , . . .)( -J-P^fe-,^, \M$( Pirs , ...)), (4.4) 

where Sj rs = (pj + p r + p s ) 2 and |A4m (pirsj • • •)) corresponds to the m-parton matrix 
element that is obtained from the (m + 2)-parton matrix element by replacing the three 
partons i,r,s by the single parton irs*. 



In general the kernel Pfff r f 3 ({zj, Sji, kj±}; e) in Eq. (|4.4D contains spin correlations of 
the parent parton. However, in the case of splitting processes that involve a fermion as a 
parent parton, the spin correlations are absent, therefore, we can write the corresponding 
spin-dependent splitting function in terms of its average over the polarizations of the parent 
fermion, 

( 8 \PqfTf.(.{ z 3> 8 jh k 3-L}'> E )W) = S ss'Pqf r f a ({Zj,Sjl};e) , (4.5) 



'According to the convention of Eq. (3.9), we have Altarelli-Parisi kernels for all the permutations of the 
flavour indices, which are obtained by permuting suitably the respective momenta in their arguments; e.g. 
Pg q g is obtained from P qgg by swapping the first two momenta. 
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where we used our generic notation for the spin-averaged splitting function (cf. Eq. (|3.8| )). 
The spin-averaged splitting function is presented in Refs. 37] explicitly. In order to bet- 
ter exhibit the soft structure of these kernels, we give different forms of the same functions. 
For the q — > qq'q' process we write 



2 ^lilrl's i Z s , S{ r , Sj s , S rs , e) 



C F T R < 



1 



Sirs S r 



z, 



+ 



Z r -\- Z s S rs (yZ r -\- Z s ^j S rs Si rs S rs (_Z r -\- Zg)^ 



Z T Z s 



Zr ~\~ Zg 



+ 



l-e 



Zr ~t~ Zg 



s rs 

Sirs 



+ (r <-> s) 



(4.6) 



We decompose the splitting functions that involve one, or two gluons into abelian 
and non-abelian contributions. For the q — > qgg splitting we only need the spin-averaged 
functions, 



P 



qgg 



c F pW + c A p qgg 



(nab) 



(4.7) 



where for the abelian part we write 



2 ^Qigrgs ( Z ii z r-> z si s iri s isi s rsj e ) 
^irs 

= C F 



1 Za 1 Ziy Zi 

+ + 



Sirs Sir SirgSig Sir Sic J %r V 1 Z- 



ez s - e(l + e)— - 



+ 



l-e 



£ - + £ )( 3 _ Zr - 2Zs ) - ^(1 - e) 

Sir Sir 



+ (r«->sU, (4.8) 



and for the non-abelian part we use 

_}_p(nab) / \ 



C F < 



SirsSr 



(l-e) 



Sir Si 



Z r z s 



Sir^s S{gZ T 
&rs (Z r + Z s ) SirsSrs Srs \%r -\- Z s )' 



-Zi 



1 



Zr ~\~ Zg 



1 (1 - Z r ) 2 + (1-Z S 
SirsSir 2 l Z r [Z r ~\~ Z s ) 



1 Zi 1 + zj 



+ 



SirSig ^2iZr Zr ~\~ Zg ^SirSr 



(l-e) 



+ 



1 



1 



2 Sirs Si 



SirsSrs 

i + (i - z, 



i + ^ 2 , i + (i-^) 2 

1 \ (z r - z s ) 



+ 



(l-e) -- 

\ %T %T 

4 2 Zg j 2 ^ 



+ 



£j- H - 



+ 



_J_( (1 _ Zl) (5: + fi_ £ 

SirsSir \ \ 



^(1 - 2 S ) \ Z. 



( £3* "T" 



+ 



SirSig 
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z r 1 + e 



[Zr + Z s 



S{ r S rs 



+ 



Zy ~\- Zg 



+ (r <-> s) > . (4. 



In the case of the g — > gqq process we decompose the spin-dependent splitting function 



as 



{lx\P gqq \v) = C F MPjffllv) + C A (fx\P^°'W) 
The explicit form of the abelian term is 



A(nab) I 



(4.10) 



1 



2 ^\Pgiq r q a i^iri Sjsj S rs , k r _\_, k s ±) \ f) 



2 

+ ■ 



1 1 



+ 



2 -Sirs 

+ K±K± " (1 " i + (r <- s) \, (4.11) 



and that of the non-abelian term is 



{l-^\Pgiq r qs (^i i%>r •,%&■> ^irt $isi $rs i &i_L ) ^V_L i ks± ) |^) 



2k- 



9' 



I w 



+ 



1 1 



SiyZ s -\- S{ s Z r Si r Si s ^ ^irs Z r Z s 



Z T -\- z s s rs [z r -\- Zg'j S rs S{ rs s rs {z T -\- Zg}* 



+ 



Si rs S r s %i 



A 7 _ hf 1 UV 
-g^ v _|_ ^' 6r ' 6s "Vs-L'VsJ 



z r -j- z s k 



9 



fJ,V 



1 



.9 



I — : r / 1 1 \ 1 1 — e z r 1 



1 



ISirgSis \ Z% Z r -\- Z s J ISiySis ^^frs ^^ir^rs H~ 

4z, 2 



+ 



+ 



2(1 - s)k^ ± K 

( Zc Zi 



7 ~~ \^s±^sA. 
Z{\Z r ~r Z s j 

+ i-e) [e rl _K ± + e sl _K ± ) 



Zi{z T -\- Zg) 



(4.12) 



where we have introduced the abbreviation 



K rs± — y „ ■ 



(4.13) 



The pure-gluon splitting kernel as given in Eq. (66) of Ref. 37] is completely symmetric 
in its indices. Therefore, without loss of definiteness, we can refer to that formula and 
obtain our definition from that by making the change of indices 1 — > i, 2 — > r and 3 — > s 
and then the formal substitution kj± — > kj±. 
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The triply-collinear factorization formulae of this section are valid in the strict collinear 
limits. If these expressions are to be used as subtraction terms, the momentum fractions Zj 
(J = i, r, s) have to be defined over the extended phase space region where the subtraction 
is defined. This definition has to be such that Zj vanishes only if the corresponding parton 
j becomes soft. Otherwise, in the kernels P qg g, Pg q q and P gg g we introduce spurious poles 
in those terms that contain 1/zj factors. The energy fractions zj = Ej/(Ei + E r + E s ) 
fulfill such a requirement over the whole phase space. Other definitions, for instance, in 
terms of invariants, like Zj = pj ■ n/[(pi +p r +p s ) ■ n], with n being a light-like momentum, 
e.g. the momentum of a hard parton k, lead to spurious poles unless the region where the 
momentum pj is collinear to the momentum n is excluded, e.g. if the subtraction term is 
defined there to be zero. 



4.2 Emission of two collinear pairs of partons 

The double collinear limit occurs when two pairs of collinear partons are emitted. The 
limit is precisely defined by the usual Sudakov parametrization Eq. ( |3.ip for the two pairs 
separately. We consider two pairs of final-state partons i, r and j, s that can be produced 
in the splitting processes ir — ► i + r and js — > j + s. The partons i and r have momenta Pi, 



p r , their parametrization is given in Eq. (3JJ. The parametrization of the second splitting 



is completely analogous in terms of momenta pj s , k^ ± , n^ s and parameter zj. 

The doubly-collinear region is identified by performing the uniform rescaling 

k r ± -> \k r ± , k s ± -> \k s ± , (4.14) 

and studying the limit A — > 0. Neglecting terms that are less singular than 1/A , the 
squared matrix element |.M m +2| 2 behaves as 



Cir-js\M^ 2 {pi,Pr,Pj,Ps, ■ ■ Ol 2 = (8vra s /i 2£ ) 2 X (4.15) 

x ( M$ (Pir ,Pjs,---)\ Pf} r (zi,z r , k rl _ ; e) Pf^ fs (zj , z s , k s r, e) \M$ (p ir , p js ,•••)) , 

where \A4m {jPiriPjsi ■ ■ •)) corresponds to the m-parton matrix element that is obtained 
from the (m + 2)-parton matrix element by replacing the two partons i and r by the single 
parton ir and the other two partons j, s by the single parton js. The c?-dimensional 
Altarelli-Parisi splitting functions were given in Eqs. Q-Q. We see from Eq. fl£lE| ) 
that the doubly-collinear limit factorizes completely as 

Cir;js\M^ +2 {pi,p r ,Pj,Ps, ■ ■ ■)? = C ir C js\M^ +2 (pi,p r ,Pj,Ps, ■ ■ ■)? ■ (4.16) 
4.3 Soft-collinear limit 

Next we consider the tree-level matrix element of m + 2 partons in the limit where the 
momentum p s of a final-state gluon becomes soft and the momenta pi and p r of two other 



-14- 



final-state partons become collinear. The soft limit is defined by the rescaling p s = X s q s 
with letting X s — > 0. The collinear limit is defined by the usual Sudakov parametrization 
of Eq. (|3,1|). The limits are taken uniformly, A r = A s = A. In this limit, the squared matrix 
element |.M m +2| 2 behaves as [37] 



CSir;s\M^ +2 {pi,p r ,Ps, ■ ■ OP = 

EE^(*) x ( 4 - 17 ) 



2 



x (MWfar, . . .)| T r Tt P<® r (zt, z r , k r±] e) \Mt ] (p ir , . . .)) 

where \-M-m ipiri ■ ■ ■)) corresponds to the m-parton matrix element that is obtained from 
the (m + 2)-parton matrix element by omitting the gluon s and replacing the two partons 
i, r by the single parton ir. In the terms of Eq. ( 4.17j ), when either j or I is equal to 
the colour charge operator is Tj r = Tj + T r and 



[ir 



2Sj(j r ) = 2(sjj + Sj r j 
SjsSs(ir) Sjs\Ssi > ^sr) 



s J(ir) (s) = —^l- = \z \ ■ i 

In the strict collinear limit 



Sj(ir){s) = Sji(s) = S jr (s) , (4.19) 

which we shall use later. Of course, CSi r - s \Ai^ +2 (pi, p r , Ps, ■ ■ -)\ 2 = if parton s is a quark. 
4.4 Emission of a soft qq-p&ir 

In the region of soft g^-emission the momenta of the quark p r and that of the antiquark p s 
tends to zero such that p r = Xq r and p s = \q s with A — > for fixed q r and q s . In this limit, 
the squared matrix element diverges as 1/A 4 and its most divergent part can be computed 



in terms of the eikonal current of a soft gluon of momentum p r + p s [37] 



Sif\M m+2 (0)( Pr ,Ps,...)\ 2 = (4.20) 

= (8na s ^f±Y, E { ^ + ^ ~ SikSrS ~ ^-P 1 ) T &\ M ™-(i *)(• • OP > 
8 " vfr,8kfr,T,s \ S i{rs)Sk(rs) S i(r S ) J 

where the matrix element on the right hand side is obtained from the original (m+2)-parton 
matrix element by omitting the quark-antiquark pair, having momenta p r and p s . 



4.5 Emission of two soft gluons 

We consider the tree-level matrix element of m + 2 final state partons when a gluon of 
momentum p r and another gluon of momentum p s become simultaneously soft. The limit 
is precisely defined by rescaling the soft-gluon momenta by an overall factor A, p r = Xq r 
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and p s = \q s , and then performing the limit A — > for fixed q r and q s . In this limit, the 
matrix element diverges as 1/A 2 and its most divergent part can be computed in terms of 
the two-gluon soft current J tllfl2 (q r ,q s ), given explicitly in Ref. |j37f| . 

The singular behaviour of | -A/i ^^_ 2 1 2 & t 0(1/A 4 ) can be written as a sum of an abelian 
and a non-abelian term [37|, 

S^\M { l 2 (Pr,Ps, ...)\ 2 = ( Sif > + S^) \M { V(Pr,Ps, • • OP • 



(4.21) 

The abelian term contains the product of two one-gluon soft functions, Sik{r) given in 
Eq. ( fUg ), 

1 m 

S^\M^ +2 (p r ,p s ,...)\ 2 = (87ra s ^) 2 - S ik(r)S^)\Mt] mm (---)\ 2 - (4-22) 

i ,j,k,l=l 

The non-abelian contribution is written in terms of Si k (r,s), the two-gluon soft function 
introduced in Ref. [37|, 

1 m 

St h) \M^ +2 (p r ,p s , . . .)| 2 = -(8vra s/ ^) 2 -C A £ S ^ s ^ \ M m,m(- • Ol 2 • (4-23) 



i,k=l 



The colour-correlated functions \^^(ik)\ 2 anc ^ 1-^ 



2 are given in terms of m- 



m,(i,k){j,l)\ 

parton amplitudes obtained from the original (m + 2)-parton amplitude by omitting the 
gluons r and s (for the definitions see Eqs. ([2.5D and (|2.6D, respectively). 



It is useful to break the sum in Eq. ( f4.22| ) into sums according to the number of hard 
partons in the eikonal factors, 

m 

£ S lk {p r )S 3l {p s )\^ 



i,j,k,l=l 



EEE E S lk (r)S 3l (s){Mt ) \{T t -T k ,T y T l }\Mt ) ) + 

i jjti k^i,j l^i,j,k 

+ 4 EE E S tk (r)S Jl (s)(M^\{T l .T k ,T j .T i }\M^) + 

i j^i k^i,j 

+ 2 EE^( r )^( s )^-i{ T '- T ^ T ^ T ^i- M ™)- 



(4.24) 



i kj^i 

The sum in the first line requires at least four hard partons, the next line involves at least 
three hard partons and the last sum contains at least two hard partons. Therefore, for the 
case of two hard partons, relevant to two-jet production in electron-positron annihilation, 
only the last sum contributes. 

For later convenience, we also record a new form of the two-gluon soft function. The 
original expression of Ref. [37| in our notation reads 



,-. / \ c.ts.0.) / \ . . Si r S ks + Si s S kr 

<S ik {r, s) = S> k '{r,s)+4: 

^i(rs)^k(rs) 



—J2 o^ik l r > S > 



—S lk (rs) , (4.25) 



- 16 - 



where 

s£°°(r, s) = S ik (s) (S is (r) + S ks (r) - S ik (r)) (4.26) 
is the approximation of the soft function Sik(r,s) in the strongly-ordered approximation 



and Sik(rs) is given by Eq. fl3,12| ), 



«5jfc(rs) 



2s 



ik 



Sj(rs) "Sfc(rs) 



(4.27) 



In order to exhibit the relation of this expression to the triply-collinear splitting functions, 
we rewrite Eq. (p~2~3D using partial fractioning, 



- ^ Sik{Pr,Ps)\M^\ i:k) ( Pl , . . . Pm )\ 2 = 



i,k=l 



l 



k^i 

SirSks Skr^is 



— (! ~ £ ) T T T T 



i 



■Sir "I - ^is ~r Sfcr T Sfcs "Sis ~\~ S kl 



2 1 

+ I + 



Sir $is $kr T" \SirSkr Sir$k. 

Sik{Sir ~i~ "Sis) Sik(Si r "I - "Sis) 

Sir(sir ^ks") SirK^ir Sis $kr $ks) 

Sik{Si r ~i~ "Sis) ( Sis Si s {Sj er $ks) 



+ 



(■Sir ~\~ >Sfc s )(sir + Si s ~\~ Sfcr Sfcs) \Si r Sk s (si r + Sis) 



+ (r <-> s) 



(4.28) 



Finally we note that if the two soft partons are of different flavours, the squared matrix 
element does not have a leading singularity, 



Srs\M { J +2 {p r ,Ps, ■■■)? = 5 frfs S r s\M ( J +2 {p r ,Ps, 



(4.29) 



4.6 Matching the doubly-unresolved limits 

If we want to regularize the squared matrix elements in all doubly-unresolved regions of the 
phase space, then we have to subtract all possible collinear and soft limits. Analogously to 
Eq. ( |3,20 ), we define the subtraction term A2|A4^ 2 | 2 , that must contain the sum 



EE 

r s^r 



i^r,s 



j=£i,r,s 



\M% +2 (pi,p r ,Pj,p s , . . .)\' z . 

(4.30) 

The symmetry factors ensure that, after summation, all collinear configurations are taken 
into account only once. The various terms in the sum ( [4.30D overlap in those regions of the 
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phase space where the limits overlap. Therefore, Eq. ( 4.30|) contains double and even triple 
subtractions. We observe that at most triple overlaps may occur because the doubly- and 
triply-collinear regions do not overlap. Let us first discuss the double subtractions. 

The triply-collinear region overlaps with the soft-collinear one. The soft-collinear limit 
{z s — > 0, z r + Zi — > 1) of the triply-collinear factorization formula, Eq. (4.4), is 



CS ir;s C iTS (pi ,p r ,p 



i7ra s ^ 2£ ) 2 T 2 



$irS(ir)s Zs 



-(M^( Ptr , . . .)\P$ r (zi, z r , k rr ,e) \M$( Pir , ...)}• (4.31) 



In the triply-collinear limit Eq. ( 4.17 ) behaves as 

(j>i,Pr,Ps, ■ ■ -)l = 

= {^a sl / £ ) 2 Tl -^—^-(M^( Pir , . . .)\Pfi{z u z r , k r± ; e) \M$(p ir , ...)). (4.32) 

SirS(ir)s Zs 

The right-hand sides of Eqs. ( |4.31| ) and (4.32) differ in subleading terms in the soft-collinear 
limit (zi r = l — z s = 1+0(A)). We also record two equivalent forms of Eqs. ( [4.31| ) and Q4.32j ) 
that will be useful later. In the i\\r collinear limit Si s = ZiS^ r )s, therefore, Eq. flOI]) can 
also be written as 



CSir:sCirs\A4 rn j r 2{Pi,PriPs-i ■ ■ -)l ~~ 

^a s ^) 2 T 2 ir -l-^(M^(p ir , . . .)\P$ r (zi, zr, k r± ;e) \M$( Pir , ...)}• (4.33) 



Sir Sis Za 



Taking the triply-collinear limit of the soft-collinear formula together with Eq. ( 4.19 ), we 



get an identical expression to Eq. ( [4.33|) , 

C irs CSi r;s \M^ +2 (Pi,Pr,Ps, ■ ■ - )| 2 = 

= (8™ s/ u 2e ) 2 T 2 r ■ ■ -)\Pf"fM, z r , k r± ;e) \M$(p ir , ...))■ (4-34) 

Sir $is %s 

The fact that Eqs. ( gj3| ) and ( p~34D are identical while Eqs. ( p~3i"D and (gT3^) differ explicitly 
by sub-leading terms is due to the ambiguity in the definition of the soft-collinear limit 
beyond the leading order, represented by Eq. ( 4.19 ). 

The doubly-collinear region also overlaps with the soft-collinear one. The soft-collinear 
limit of the doubly-collinear formula, Eq. (4.15), is 



CSir;s(I!j r .y S |.A4 m ^_2(Piji'r ,Pj ,p s ,...) 



(4.35) 



Tra^fT 2 -J—±(M$(pir, Pj , . . .)\Pj?Uzi, z r , k r± ;e) \M$(p ir , Pj , ...)>, 



while in the doubly-collinear limit Eq. (4.17) behaves as 



Cir;js££*ir;s I -^m+2(^ i Pri Pj i Psi • • •) 



(4.36) 



lira s /j, 2e ) 2 T 2 



(M$(pir,Pj, . . .)\P)^ r {Zi, Zr, k r ±;e) \M^>(pir,Pj, ...)}. 



,(0) 



(0), 
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The difference between Eqs. ( 4.35| ) and ( 4.36 ) is again subleading in the soft-collinear limit 
( Zj = l-z s = l + 0{\)). 

The region of splitting into a collinear /jgg-triplet (parton i can be either a quark or 
a gluon) overlaps with the soft gg-emission. In this overlapping region, 



(87ra sM 2£ ) 2 T?r R - 



C irs Slf>\M ( Z2(Pi,Pr,Ps,---)\ 

2 



(4.37) 



%s Si(rs)Srs(z r -\- Z s )' 



In the soit-qq limit, this formula differs by a subleading term from the doubly-soft limit of 
the triply-collinear factorization formula, 



,(0) 

37ra s/ / 2£ ) 2 T 2 T R - 



Si,fC irs \M^ +2 ( Pl ,p r , Ps ,...)\' 

2 



(4.38) 



l (ySi r Zg Si S Z r ^ 

Z r -\- Z s Si( rs yS rs (z r -\- z s 



\M^( PlJ ...)\' 



The soft (^-emission does not overlap with any other doubly- unresolved regions. 

The phase space regions where soft gluon-pair emission occurs overlap with all other 
doubly- unresolved parts of the phase space (except for the soft-gg emission) . In the over- 
lapping region of splitting into a collinear /j^g-triplet (parton i can be either a quark or a 
gluon) and soft ^-emission we obtain 



C irs S^\M^ +2 \ 2 = C irs S^\Mt\ 2 \ 2 + C irs SiT b) \M { Z 2 \ 2 , (4.39) 



where 



C irs S^\M^ +2 ( Pi , Pr , Ps ,. 
C irs SiT h) \M^, 2 ( Pl ,p r ,p s . 



^a sM 2£ ) 2 T 2 T 2 . 



Azf 



Sir^is^r Zg 



■\M^( Pi ,...)\\ (4.40) 



Zi 



Si(rs)Srs Si(rs)$rs{Zr ~\~ Zg)^ Si(rs)$rs \Zf H~~ Z s Z< t 



2zf 



1 



s i(rs) s ir %r (z r + Zs) S { 

(rs) s is z r (z r + Z s ) 



+ —^— ( — + ) +{r^s) 

SirSrs \Zs Zr H~~ Z< 



\M^( Pu ...)\ 2 . 



(4.41) 



Computing the doubly-soft limits of the Pjigg splitting kernels, t we can derive that 



^Eq. (4.11) and the azimuth-dependent terms in Eq. (4.12) and in the triple-gluon splitting function do 
not have leading singular terms in the doubly-soft limit, as expected. 
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S^C irs \M^ + 2(Pi,Pr,Ps,...)\ 2 

= (^a s ^) 2 T 2 {T 2 



+ 



+C A 



1 



1 



(1 6") {Sir^s Sis^r) 

2 1 11 

+ 



(z r + Z s ) Si (z r + Z s ) 



+ - 



1 



1 

— + 



1 



+ (r <-» s) 



•|^)( Pi ,...)| 2 . (4.42) 



We see that in the doubly-soft limit Eqs. ( |4.37 ) and ( 4.38Q differ by subleading terms, 
as well as the sum of Eqs. (4.40) and ( |4.41| ) differ from Eq. ( [1.42 ) by subleading terms 
( Zi = l- z r - z s = l + 0(A)). 

In the doubly-collinear limit, the doubly-soft factorization formula simplifies to 



C ir . js S^>\M^ + 2(Pi:Pr,Pj,Ps, •••)! = 

= (8,a s ^) 2 T 2 A * T 2 2l Z ± \M^( Pl ,p 3 , . . .)| 2 



(4.43) 



In the region of the doubly-soft gluon emission, the limit of the doubly-collinear formula, 
S^C ir . JS \M^ +2 (p u p r , Pj ,p s , ...)\ 2 = 



(87ra s ^ 2e ) 2 T 2 



|M#Wi,-)l 2 , 



(4.44) 



differs from Eq. ( 4.43 ) by subleading terms (zi = 1 — z r = 1 + 0(A), Zj = 1 — z s = 1 + 0(A)). 

The last doubly-overlapping region is the overlap of the soft-collinear region with the 
emission of two soft gluons. In the soft-collinear limit Eq. ( 4.21 ) becomes 

CSi r;s S^)|^ 2 (pi,p r ,p a ,...)| 2 = 

= -(8na s ^) 2 A ^T 2 £ £ k,(-)|A<£*>' . . .)| 2 . (4.45) 



i ¥s 



Note that the non-abelian part of Eq. (4.21) does not yield a leading singularity in this 
case. In the same region, the doubly-soft limit of the soft-collinear formula is 



S^CS ir;s \M^ +2 {p i: Pr:Ps, 



2c 



(4.46) 



3 ¥3 



Thus, in the doubly-soft region, Eqs. ( [4.45| ) and fl4.46|) differ by subleading terms 
1-Zr = 1 + 0(A)). 
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In addition to the pairwise overlapping doubly- unresolved regions, there are two regions 
where three limits overlap. The first one is the overlap of the triply-collinear, the soft- 
collinear and the doubly-soft emissions. In this region, 



Cj rs CSj r .; S S rs |A'(^j_|_2(pi)Pr)i'si ■ ■ ■)! — 

= (8na s ^) 2 - — T 2 — — T\ |m£) (p i} . . .)| 2 

Sir %r $is 



(4.47) 



The other region of triple overlap is that of the doubly-collinear, the soft-collinear and the 
doubly-soft emissions, where 



Cj r; j s CSj r;s S rs |.A/i| n ^_2 (Pi'Pr)PjiPs) ■ • 01 — 

= (8na s ^) 2 - — T 2 — ^ Tf \M^{p^ . . .)| 2 

Sir %r Sjs %s 

Using Eqs. (jDfi])-(jUSi), we can now easily check the following relations, 



CS«r;sfCj rs — Cj rs CSj r;s ^ |-Mm+2 (P*'PsiPr 
CSjr;s (^^ir;js Cj r -j S CSj r ; S ^ l-^ / ^m+2(P*' Pr> Pj : i Ps 

CSjr;s (Cir-Js^rs Cj r ;j,sCSi r ; S Srs^ \^^m+2^P^^ P r ' Pj ' P s 

Srs ^Cj rs — Cj rs S rs ^ j ■M-\ r ^j r 2 {Pii Pr > Ps 

^rs(^-^ir;js Cj r -j S S rs ^ (Pi) ; Pj : ) Ps 

S rs ^CSj r;s — CSj r;s S rs ^ [ -M^^ (Pi) Pr ) Ps 

,(0) 



Srs ( Cj rs CSj r ; S Cjj-sCSjj-jsSf.^ J j-'^ ?n _|_2 (Pi ) Pr > Ps 



SrS ( C^jj'gCSjy.'; 



Cjr-jjsCSjy-jsSf-s j \-M-\ r ^j r 2^\Pi^ Pr > Pj > P 



)| 2 
)| 2 
)| 2 
)| 2 
)| 2 
)| 2 
)| 2 
)| 2 
)| 2 



(4.48) 

(4.49) 
(4.50) 
(4.51) 
(4.52) 
(4.53) 
(4.54) 
(4.55) 
(4.56) 
(4.57) 



With Eqs. ( |4.49| )-( |4.57| ) at hand it is not difficult to verify that for the subtraction term, 

r s^r L iy^r,s L j=£i,r,s 



+ 2 ^rs 



ij^r,s 



C&ir;s^rs 



j^i,r,s 
1 _ 



+ 



C£>ir:s£>rs J -)- 



x \M^ +2 (pi,p r ,Pj,p s , . . .)| 2 , 



(4.58) 
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the following equations hold, 





-A 2 )|^ +2 | 2 = 


0, 




-A 2 )|-M^ 2 I 2 = 


o, 


CSj r;s (l 


-A 2 )\M«U 2 = 




S rs (l 


-a 2 )|^ 2 | 2 = 


0. 



(4.59) 
(4.60) 
(4.61) 
(4.62) 

Therefore, Eq. ( 4.58j ) is free of double and triple subtractions. Instead of Eq. Q4.58 ), the 
subtraction terms that can be derived from antennae factorization, as proposed in Refs. [16 



29, BOB, can also be used for regularizing the infra-red divergencies in the doubly- unresolved 



regions. The advantage of antennae subtractions over our proposal is that those avoid the 
problem of defining the momentum fractions. Note however, that contrary to the original 
antennae factorization terms of Ref. [22], the spin-polarization of the parent parton in the 
collinear regions is not yet taken into account in Refs. [16, |29|, 3C[. The spin-averaged 
antennae terms can only be used in a hybrid slicing-subtraction method as applied in 
Ref. H|. 

We can simplify Eq. ( |4.58| ) by removing the terms that contain triple limits. Comparing 
Eqs. flOa|) and ( fU8|) we see that in the last term in Eq. ( |4.58|) 



J ir;js 



Is 

2 rs 



\M {0) I 2 



^ ■ Is \M {0) I 
y f ysn rs l y m+2l 



(4.63) 



Furthermore, comparing Eqs. flPcD and flOlD to Eq. (^47|), and using Eq. flQ7|), and 
observing that the soft-collinear limit of Eq. ( [4.41| ) is zero, we see that 



S rs ab) - Srs b) ) 5f r g 5f s g + S rs 5f r g 5f s g \M 



_i c . 

This equation shows that we can simplify the notation by introducing 

c(A) _ c(ab) r e 
°rs — 3 rs u f r g °f„g 

and 

= SiT b) 5 fr9 5 fsg + S rs 5 f „ a 5 



m+2\ 



Irs Of rq Of s g . 



(4.64) 

(4.65) 
(4.66) 



With this notation S r , 



5. Matching the singly- and doubly-unresolved limits 

Eqns. ( p59H462^ ) ensure that the difference \M^ +2 \ 2 — A2\Aif^ +2 \ 2 does not contain leading 
singularities (therefore, it is integrable) over those doubly-unresolved parts of the phase 
space where the factorized matrix elements \Mm ) m the subtraction terms are finite. 
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Nevertheless, the regularized squared matrix element, |.M^j_ 2 | 2 — A 2 |.M^:, 2 | 2 , contains 
leading singularities over the singly-unresolved parts of the phase space. We have a simi- 
lar problem concerning the subtraction term Ai|.M^ 2 | 2 that is used for regularizing the 
singly-unresolved limits, namely, Ai|.M^ 2 | 2 becomes singular over the doubly-unresolved 
parts of the phase space. In a general subtraction scheme one wishes to regularize these un- 
wanted singularities, introduced by the subtraction terms, by defining further subtractions 
that take into account the overlap of the singly- and doubly-unresolved regions. The only 
consistent way to achieve that is to introduce a third subtraction term — Ai 2 |.M^^ 2 | 2 that 
simultaneously cancels all unwanted singularities in both the singly- and doubly-unresolved 
parts of the phase space. As a result, the combination 



(A 1 + A 2 -A 12 )|^ 2 I 



(5.1) 



may serve as a NNLO subtraction term that regularizes the squared matrix element in all 
relevant unresolved regions of the phase space, i.e. the following equations must hold, 



C ir (Ai 
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-A l2 )\M% 


2_ 


Cir|-M m+2 | 2 , 


S r (A a 


+ A 2 


-A 12 )\M^ +2 


2_ 


S 1 A// (0) I 2 


Cj rs (Ai 


+ A 2 


-a 12 )\m£> +2 


2_ 


Cirs\-M-m+ 2 \ 2 ! 


Cir-Js ( Al 


+ A 2 


-A 12 )\M^ +2 


2_ 


C ■ \M {0) I 2 


GSir;s(Ai 


+ A 2 


-A 12 )\M^ +2 


2_ 


CS ir;s |A^ 2 | 2 


Srs(Ai 


+ A 2 


-A 12 )\M^ +2 


2_ 





(5.2) 
(5.3) 
(5.4) 
(5.5) 
(5.6) 
(5.7) 



In order to prove Eqs. we define Ai 2 |.M^ 2 | 2 as AiA 2 |.M^ 2 | 2 . Explicitly, 



M2\M^ +2 \ 2 



E 



S t A 2 + ~ C kt A 2 - ]T C kt S t A 5 



\M (0) I 2 

\ Jvl m+2\ > 



(5.8) 



where 



s,a 2 |^ 2 i 2 =j; y, 



Srv o go c( na t') 



+ S t S rt \\M^ +2 \\ (5.9) 



C fct A 2 |^ 2 | 2 = Crf + £ 



+ CktCktr + CktCSkt;r — (5.10) 



(N) 



■l^ 2 l 2 , 
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and finally 



(0) |2. 



C k & t A 2 \M^ +2 



a c (nab) 



+CjttSfCfc r t + CfaStSrt — CkfStCSkr-tSrt 



CUC c( nab ) 



(5.11) 



Note that the direct calculation of the terms on the right-hand side of Eq. (|5^) yields 
significantly larger expressions than those in Eqs. (gTg)-(gTTg) . In Sect.^ we shall derive 
some factorization formulae we used in obtaining Eqs. (|5.9|)- (|5.1l| ). In fact, all the terms 
of Eqs. (|5T9D (^7Tl|) are singly-unresolved limits of the factorization formulae for doubly- 
unresolved emission. Most of them can be computed easily using the explicit expressions 
of the previous sections. The only non-trivial ones are those that are related to the iterated 
singly-unresolved limits, which we discuss first. 



6. Iterated singly-unresolved limits 



In this section we derive the factorization formulae valid in the singly-unresolved limits 
of Eqs. (gig) and ( gl^) that we call the iterated singly-unresolved limits. Thus we give 
the explicit expressions of the various terms in Eqs. (|3.23p and ( 3.24 ). All the factorization 
formulae presented in this section are valid for helicity- and colour- correlated squared matrix 
elements in Eqs. (3. It) and (3.1$) written in terms of on-shell momenta. 



6.1 Collinear limit of the soft-factorization formula 



In the derivation of the collinear limit of the colour-correlated squared matrix element 



l^ro+i-fj I) I 2 °^ ( |3 . IIP , we have to distinguish two cases: (i) the partons of the collinear 



pair i, r are different from partons j and I; (ii) one parton of the collinear pair is either j 
or I, the other is different. In principle there is a third case, when j and I become collinear. 
The corresponding 1/s,-; pole is suppressed by the sji factor in the eikonal factor that 
multiplies l-M^+i-Q' z)l 2 > t nereiore ; the contribution is subleading. Also, [A^^./- « | 2 does 
not have a leading singularity when i and r are both quarks or antiquarks. 



In the first case, the colour-charge operators do not influence the computation of the 
collinear limit, l-M^^,. «| 2 factorizes on the collinear poles as I-M^+J 2 does. In the 
second case, due to coherent soft-gluon emission from unresolved partons, only the sum of 
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I^Cl;W)| 2 + \M { X,(r,l)\ 2 ( OT \ M Zv, U J 2 + M£Uy,r)| a ) faCt ™ 



as 



S ii( s )l^ 1;W )| 2 +^WN!li ;W) | 2 * 8vra s ^5 (4r) K S ) (A^l T ir ■ T, \M$) , 

(6.1) 

where TV = Tj + TV and Sr ir \j(s) = Sj( ir }(s) is denned in Eq. ( [4.18 ). Collecting all the 
terms, in the collinear limit i\\r we obtain 

C ir S s \M^ +2 ( Pt , Pr , Ps ,...)\ 2 = (6.2) 
= -(8vra s /i 2e ) 2 x 

x 7-| E ^ jl (s)(M^\p ir ,...)\Pi%Xz l ,z r )T r T l \M^(p ir ,...)) + 
+ ^ 5 (ir) K«)(A^ ) (Pir,..0l^3r(^^) T ^- T H^ ) (Piri--0>}, 



that is, 



#(ir) 

C 4r S s |^+ 2 | 2 = CS-; S |^ 2 | 2 . (6.3) 
6.2 Soft limit of the soft-factorization formula 



We are interested in the limit of the right hand side of Eq. (3.11) as a gluon of momentum 



p s = Xq s becomes soft (A — > at fixed q s ). The soft limit of the colour-correlated squared 
matrix elements, fc) I 2 ' can De obtained in a similar way as that of the squared 

matrix element. We have to distinguish two cases: (i) s is different from i and k, (ii) either 
i or k gets soft. 

Let us start with case (i), when we can immediately use the well known factorization 
of the matrix element \A4^ +1 (p s , ...)}, 

(u\M^ +1 (p s , ...)>=< g B ifJ v { Ps )\Mg){. . .)} , (6.4) 

where we neglected all contributions that are less singular than 1/A. The factor J u (p s ) is 
the eikonal current for the emission of the soft gluon s, 

r{ Ps ) = Y^T 3 ^. (6.5) 

We substitute Eq. ( |6^ ) into \M^ +l .^ k ^{p s , ■ ■ -)\ 2 and obtain 

^\M^ +1 . m ( Ps , . . .)| 2 = -8vra s ^ \SM<M&K ■ -)\TfTtT b k Tt\Mt ] {. ■ •)) • (6.6) 

3,1 

We can use the colour algebra to derive the identity 
2T a T b T b T a = {T . . TuT . . Tk} _ 

: £ / rri Cirri Crrib £ i rfiCLrpbrpC C i rjiarjiCrrib C i rjia rnb rriC £ \ i 

-Uabc{lj J-i +k°li + J-j J-i J-k^lk + L l J-i IkOjk + J 1 J-i J-k°ji) + 

+CATiT k (5u5j k - 5i k 5 jk - 5u6ji + Si k 5ji) , (6.7) 
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which we substitute into Eq. fl6.6j). Performing the summation over j and I, we find that the 
terms proportional to f a b c in Eq. ( |6.7[ ) are symmetric in b and c, while f a b c is antisymmetric, 
therefore, the sum of these terms is zero. The remaining terms yield 



^\M^ +1 ; { i,k)(Ps,---)\' 



-87TQ S /X 2£ - 



3,1 



(0) , 
m;(i,k)(j,l) ' 



+ 2C A S ife ( S )|A^ ( . fc) (. 



.(6.8) 



In the second case, when for instance gluon i gets soft, we first use colour conservation, 



,(o) 



(6.9) 



We obtain the factorization formulae for the \M m+1 .y k ^ | colour-correlated matrix elements 
the same way as we did in case (i), and substitute those into Eq. {J6.9|) , 



.10) 



In the first term we can perform the summation over I using colour conservation, and find 
that the sum is zero. Therefore, 



S.|A*S. l!(M ,(P- 



8na s ^lc A Y,Si k (s)\M^ m (...)\ 2 . (6.11) 



We check explicitly the validity of Eqs. ( |6.8j ) and ( |6.11| ) in Appendix A using the known 
colour-correlated squared matrix elements for the e + e~ qq + ng (n = 1, 2) processes. 

We can now easily compute the soft-factorization formula for the right-hand side of 
Eq. 03.1l|) . We decompose it as 

j,fc i,k^s k 

In the soft limit all the eikonal factors are finite, so we can simply substitute Eqs. (| 
and (|6.11|), and obtain 



S s S r \M^ +2 ( Pr , Ps ,...)\ 

,1 



(8Tra s /j, 



2e\2 - 



E ^*(r)5,-«W|A4^ )M (. 



i,k,j,l 



i,fc=i 



(6.13) 
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where 

s£°'\r,a) = S ik (s)(s is (r) + S ks (r) - S ik (rj) (6.14) 

is the two-gluon soft function in the strongly-ordered approximation. If either r or s is a 
quark, then the result is zero. 

Finally, we consider the region where the soft momentum p s becomes also collinear to 
the momentum of another parton j. In this region, we have 

ls jk (s)ls js (r)+S ks (r)-S jk (r)) * Sjk ^ 0(1/X 2 ) , (6.15) 

thus the term proportional to Ca is subleading. The eikonal factor Sji(s) becomes inde- 
pendent of I, 

M s ) - — - > ( 6 - 16 ) 

therefore, using colour conservation, we can perform the summation over I in the first sum 
on the right-hand side of Eq. ( |6.13 ), and obtain 



C js S s S r \M^ +2 ( Pr ,p s , ...)\ 2 = ~{^a s ^f^ Z fT) t;Sik(r)\M% (. . .)f . ( 6 .17) 

Comparing this result to Eq. (|6.2D , using Eq. ([T^) and 

S s CS js . r \M^ +2 ( Pj ,p r ,Ps, • • -)| 2 = 8ira s ^ — - T) S r \M^ +1 ( Pi ,p r , . . .)| 2 , (6.18) 
we see that 

S s (c js - C JS S s )s r \M^ +2 (p r ,p s , . . .)| 2 = 0. (6.19) 
6.3 Collinear limit of the collinear-factorization formula 



Next, we consider the limit of Eq. ( 3.10| ) when parton s becomes collinear with another 



parton j. The collinear-factorization formula involves the helicity-dependent tensors, 

T$i = (M^ +1 ( Pir , Pj , Ps ,..m(h'\M { Zi(Pir,Pj,Ps, •••)>, (6.20) 

which are the leading-order (m + l)-parton squared matrix elements not summed over the 

w - \ M (°) 12 



spin polarizations of parton (ir), i.e., if parton (ir) is a gluon, then —g^Tj^^ = |.m£ 2 



and 5 S s'T~m+i = l-^m+il 2 ^ (* r ) i s a q uar k- 

The computation of the collinear limit of the T^f^ tensors follows that for the collinear 
factorization of the squared matrix element. When j ^ (ir), we obtain 

C js %!^l(Pj,Ps,Pir,---) = (6.21) 

= 8™ s/ u 2e — (M$(p js ,Pir, ■ ■ .)\h) Pf\{zj,z a ,k a r,e) {h'\M$( Pjs ,Pir, ...)>■ 
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If j = (ir), then we find 
C(ir)s^m+l(Pir>Psi ■ ■ ■) 



(6.22) 



= 8na s fi 2e ^—(M$(p {ir ) s , . . .)\ P^' fs (z s , k s r,p {ir)s ,n {ir)s ; e) \M ( ^(p {ir)s , ...)}• 

The Pffi(z,kx ;p,n;e) kernels (splitting tensors) are related to the ci-dimensional Altarelli 
Parisi splitting functions by summation over the spin indices h and h', 

jhh 



$hh> {s\P q h 9 h ' (z, k ± -,P, n; e)\s') = (s\P$(l - z, z; e)\s') , 
S h h'(^\Pqq h '(z,k±;p,n;e)\iy}= (n\P$(l - z,z,k±;e)\v) , 
d a p(jp,ri){s\P]]£ (z,k±;p,n;e)\s'} = (s\P$(z, 1 - z;e)\s') , 
d a p(p,n)(fJ,\P£f(z, 1 - z,k±;p,n;e)\u) = (n\P$(z, 1 - z, k±;e)\u) . 

Here d a p(p,n) is the gluon polarization tensor for physical polarizations, 

Pan/3 + Ppn 

a 



d a p(p,n) = -g af3 + 



p ■ n 



(6.23) 
(6.24) 
(6.25) 
(6.26) 

(6.27) 



In order to obtain the collinear limit of the collinear-factorization formula, we insert 
Eqs. dOID and ([[22]) into Eq. flgip . For j ^ (ir) we obtain 



C js C ir \M^ +2 (pi,Pr,Pj,Ps, • • OP 



57ra s /i 



2e\2 



(6.28) 



X (M$ (p ir , Pjs, ■ ■ •) I P)°) p fcp± ! e) Pfff s (z j ,Z s ,k j ±-,£)\M$ (p ir , Pjs ,•••)) ■ 

The quark splitting function, (/i|P<^ z g , k±; e)\h'), is diagonal in the spin indices (see 
Eq. (jO|)), therefore, 

{h\Pj$(z q ,Zg,k r ±;£)\ti)P^(z s , k s± ;p {ir)s ,n {ir)s ;e) = pW (z q , %;e)P g ( J](l - z s ,z s ;e) . 

(6.29) 

As a result, the factorization formula in the case of parton j = (ir) being a quark is 
C^ r)s C ir \M^ +2 (pi,p r ,p s , • • -)| 2 = 



) PJ? (zg , z 9 , /c r j_ ; e) x 



S 7TQJsM ) ' </./ 1 -'/' '//• 

s ir s (ir)s 



x(M$(p (ir)s , . . .)! P$(l - z s , k ir ±;e) \M$(p (ir)s , . . .)). 



(6.30) 



Thus we need explicit expressions only for the P^? splitting tensors, 
(s\P^(z,k ± ;p,n;e)\s') = 



C F 5 S 



— ff^-2 Z kl ^ + (l- z ] 



1 - z k\ 



pa n (3 _|_ pP n ot 
1p • n 



+ 



(^\Pgg(z,k ± ;p,n-e)\v) 



2C A 



1 - z 



•far" + 



A/ I fb I 



1-Z 



z(l - z)d a ?(p, n) 



fV I (V I 



+ 



(6.31) 



(6.32) 
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In Eqs. ( |6.31| ) and ( 6.32 ) the ellipses mean neglected terms that do not contribute to the 
contraction of the splitting functions and splitting tensors 



(a | P$ (zi,z r ,k r± ;e)\P) P^ (z s , k s ± ; p {ir)s , n {ir)s ; e) 



(6.33) 



(summation over the gluon polarization indices a and f3 is understood), which appear in 
the factorization formula for parton j = (ir) being a gluon, 



J (ir)s 



Cir\M% +2 (Pi,Pr,Ps, ■ ■ -)\' 



57ra s /x 



2e\2 



Sir 8(ir)s 



— (a\Pjig>{zi,z r ,kr±-,e)\P) x 



X {M$ (p {ir )s >---)\ P gf s ( Z s, h± ; P{ ir)s ? *^{ir)s ' ^ 

)\M$(p {ir)s ,...)). (6.34) 

In Appendix B, we use the known T^ v tensor for the process e + e~ — > qqg to check the 
validity of Eqs. (|6.3l| ) and ( 6.34 ) explicitly. In Appendix C, we translate Eq. ( 6.34| ) to a 
helicity basis and compare the obtained expression to the corresponding one derived using 
the known collinear factorization properties of helicity amplitudes. 



We note that the collinear-factorization formula for the tensor in Eq. ( 6.22j ) fulfills 
gauge invariance (only) in the collinear limit, 

( <Mg?(p(ir)« • • -)\ P 9fS z > k ^P,n-e)\M^{ P{ir)s , ...)) = 0(fcj_) , 



Pi 



3.35) 



where pf r = (1 — z)p? ir ^ s — kj_ + 0(k^_), and in obtaining Eq. ( |6.35| ) for f s = g, we used 
P{ir)si^m = 0. We kept the gauge terms in Eqs. ( |OlD and ( |6.32| ) only to maintain 



this gauge invariance; those terms do not contribute to either Eqs. ( |6.25 ) and ( |6.26| ) or 
Eq. O. 



6.4 Soft limit of the collinear-factorization formula 



As in the case of the soft limit of the soft-factorization formula in Sect. 6.2, in order to 



obtain the soft limit of the collinear-factorization formula we use the factorization of the 
matrix element, Eq. ( |6.4| ), which we substitute into Eq. (3.4). If the parent parton of the 
collinear splitting is a gluon, the factorization formula contains non-trivial spin correlations 
(see Eqs. (^) and fl3.£Q ), therefore, we have to distinguish two cases: (i) (ir) is a hard 
parton; (ii) (ir) is the soft gluon. If the soft gluon is labelled with s, we find that 



if (ir) is a hard parton, then 

S s Ci r \M^ +2 {Pi,Pr,Ps, ■ • OP 
= -(8vTQ s /i 2£ ) 2 X 



+ S {ir)l (s)(M<$(. 



p^Tj.t^mM (...)) + 



Pf^-T^M^ (...)) ; (6.36) 
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if (ir) = s is the soft gluon, then 
S(i r )Ci r \M^ +2 {pi,Pr,---)\ 2 = 



(6.37) 



where 



so g^ v S^(s) 
Eq. «, 



qfJ,Ur 



(6.38) 



Sji(s). Notice that the right-hand side of Eq. (|6.36| ) coincides with that of 



S s C ir \M^ +2 \ 2 = CS ir -s\M^ +2 \ 2 



(6.39) 



If the soft momentum p s becomes collinear to the momentum of another parton j, then 
Eq. ( |6.36 ) becomes 



Cj s S s Ci r \M m+2 (pi,p r ,p s ,Pj, ■ ■ ■) 



57ra s /i 



2e\2 



1 2 



i2 *3 



{M^{...)\py> fr \M^ {■■■)), 



(6.40) 



where j can also be (ir). 



If (ir) = s is the soft gluon, then in the j\\s limit Eq. ( 6.37 ) becomes 

Cj(i r )S(j r )Cj r |.A/f^|_ 2 (pi,p r , . . .)| 2 = 

= (8vra s/ ^) 2 - — T ] — 2C\ (— + — ^ \M$(. ■ 
if i and r are both gluons, and 

= (87ra s /x 2 



(6.41) 



■<<->- ' 2 T 2 ^Ir|MW(...)| 2 



■Sir ^j(ir) ^ir 



(6.42) 



if i and r are quark and antiquark, where we used k± ■ p^ = 0. Comparing Eq. ( S.28| ) to 
Eq. ( |6T4Q[) and Eq. (|oT3l ) to Eqs. (16.41|) and ful ), we find that 



S s (c is - C js S s )c ir \M^ +2 (pi,p r ,Ps,Pj, ■ ■ -)| 2 = 0, (6.43) 

no matter whether or not the soft parton s is (ir), or whether the soft parton becomes 
collinear with the hard parton (ir). 

We remark that the treatment of colour in the derivations of these iterated singly- 
unresolved limits was different for the spin-correlated matrix elements in Eq. (|3.4j ) from that 
for the colour-correlated matrix elements in Eq. ( |3.11 ). This means that soft factorization 
formulae do not exist for the simultaneously spin- and colour-correlated squared matrix 
elements. 
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6.5 Doubly-unresolved limits in the strongly-ordered approximation 



The factorization formulae of the iterated limits coincide with those valid in the strongly- 
ordered doubly-unresolved parts of the phase space. This statement can be used in proving 



Eqs. (5.2)-(pJ), therefore, we discuss the strongly-ordered limits briefly. 



First we consider the strongly-ordered doubly-collinear limit, which can occur in two 
ways. One is defined by the Sudakov parametrization of the momenta of two separate 



splitting pairs ir — > i + r and js — > j + s as in Sect. 4.2, but with the strongly-ordered 
rescaling, 

k r ± -> X r k r ± , fc s _L -> X s k sl _ , (6.44) 

and studying the limit A r , X s — > 0, with A s 3> A r . Keeping terms that are of 0[1/(A^A^)], 
the squared matrix element |7W m+ 2| 2 behaves precisely as described by Eq. (|4.15| ), 

The other strongly-ordered doubly-collinear limit occurs when there are two subsequent 
collinear splittings from the same hard parton, irs — > ir+s — ► i+r+s, which is the strongly- 



ordered limit of the triply-collinear emission discussed in Sect. |4.1| . In order to describe the 
ir — > i + r collinear splitting, we introduce the following Sudakov parametrization of the 
parton momenta, 

A = (l-z,)4 (6.45) 

where pf? r is a light-like momentum that points towards the collinear direction of partons i 
and r, nf r is an auxiliary light-like vector that is orthogonal to the transverse component 
of the third parton (rii r k s ± = 0) and n are the momentum components that are orthogonal 
to both pi r and rii r . Without loss of generality, we may choose nj r = nu r \ s . The other 
vectors are defined in the same way as after Eq. ( |4.1| ) . Comparing the two parametrizations 
of momenta p, and p r in Eqs. Q and ( ggg ), we deduce that 



Q = -^, j = i, r (6.46) 

and 

kix = -Qk s ± - k , k r x = -( r k s ± + k . (6-47) 

The strongly-ordered collinear region is identified by performing the rescaling 

k — ► Xk , k s ± — > X s k s ± (6.48) 

and studying the limit A, X s —> 0, with A s 3> A. In this limit Eqs. (|4.4j )- fl4.12|) , as well as 
the triple-gluon splitting kernel take much simpler forms. In fact, when neglecting terms 
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that are less singular than 1/(A 2 A 2 ) the squared matrix element |.A/f m +2| 2 behaves as 



C( ir ) s C ir |A^^ 2 (pi,p r ,p. 



57ra s /i 



2e\2 



1 1 



Sir S 



(ir)s 



(6.49) 



X {pirs-, • • •) I Pfifrfs CC* j Cr: K ) 2 s> &S-L) s ir) Sis? s rs! s) (Pirsj • • •)) 



where the kernels Pf°f r f a are those in Sect. 4.1 in the strongly-ordered approximation. 
Explicitly, 



Pqsq'qL (Cii Cr , Z s , Sir, Sis, S r s', 



R 



P(°)(z s , 1 - e) - 2C F Q ( r (1 - z s ) + 



"^9iSrSs (Ci I Cr j | £) = -P<jg s (1 — z s, z s\ e )Pqig r (Ci ! Cr i £ ) ) 
Pq s gig r (Ci; Cr; K ; 2 S , Sir; s is; s rs; e ) = 



(6.50) 

(6.51) 
(6.52) 



2C A 



P^fe, 1 - z s ;e) + C F (1 - e) (& Cr (1 - z.) + 

\Sr y / V s lr S(i r ^ s J 



where s KS = 2k ■ p s . Furthermore, 

^l^r(Ci,Cr,^,fc s± ;e)|^) = P^CCr.CijeXMl^?^-,! - , (6-53) 

(lAPg s qiq~ T (Cii Cr, K , Z s , k s ±, Si r , Si s , S TS ', E)\l') = 



2CaTr 



1 - 



+ sT — — + 4&Cr 

SirS^ir^s z s K 



-4C A (1 - e) Zs (1 - 2.)Pii(Ci, Cr, «; e) ^ - 



(6.54) 



(A'l^giSrSs (Ci) Crj K ) ^s_L) s ir j s is) s rsj 

2.9 1 - 2. 



4C 2 . 



9' 



1 



Ci , Cr 
Cr Ci 



1 — £ H 2 

-AIV KS 

9 



1-2, K^/C" 



-2(1 - e) Q £ 

Z s K 

-4C7 A (1 - e) z s (1 - Z S ) P<f (Ci, Cr) • 
In actual computations we can use 

«" = CrPf " + 0(K 2 ) , fc^ ± = Z sV l - Z irP » + 0(k 2 ± ) . 

The Pf °t t kernels can also be written in a unified form, 

Jijr Js 

(X^fifrfs (Cii Cr , z s, k s ±, Si r , Si s , S rs ', £)\x) = 



(6.55) 



(6.56) 



- 32 - 



= ^^I^A/^^' ^' ^' ^I^') <XI^^>.(^' ^^i^,^^; eXxO , (6.57) 
h,h' 

where we used (see Appendix B) 



40 Cr *' {K ' 2 % ±)2 = + 0(k s± ) . (6.58) 

1 Z s K SirS(ir)s 



Therefore, Eq. coincides with Eqs. ( fOg ), ( ]6~3C| ) and ( |B~3^ ). 

The strongly-ordered doubly-soft emission is defined by parametrizing the momenta 
p r and p s as p r = X r q r and p s = X s q s , and taking the limit A r , X s — > 0, with X s ^> X r . 
In this limit the squared matrix element |A4 m+ 2| 2 does not have a leading 0[1/(A^A 2 )] 
singularity when a soft gg-pair is emitted. In the case of the emission of two soft gluons, 
the singular behaviour of |.M m +2p can be written as in Eq. (B) , the sum of an abelian 
and a non-abelian contribution, 

1 / 1 m 

S s S r \M^ +2 ( Pr , Ps , ...)\ 2 = (8na s ^) 2 4 U E ^(^) 5 ^fe)I^S^)0-,ol 2 " 

V i,j,k,l=l 

m \ 

-CaI:^ 00 ^,^)^^)! 2 , (6-59) 

i,fc=l / 

with defined in Eq. ( gjg ). Thus we see that the abelian term is the same as in 

Eq. ( ^24] ), while the non-abelian term in the strong-ordered limit is simplified a lot. We 



also see that the strongly-ordered doubly-soft limit is the same as the iterated doubly-soft 
limit recorded in Eq. ( |6. 13 ) . 



We remark that all the formulae in the strongly-ordered doubly-soft limit are symmetric 
in r and s, therefore, in the limit A r ^> A s , exactly the same formulae are valid. 

Finally, we discuss the strongly-ordered soft-collinear limit, which is defined using a 



similar parametrization of the momenta as in Sect. 4.3. The soft limit is defined by the 
rescaling p s = X s q s with letting A s — > 0. The collinear limit is defined by the usual Sudakov 
parametrization of Eq. (|Q|), with rescaling of k r ± as k r .± = X r k r .± and letting A r — > 0. The 
limits X r , X s — > are taken such that either X s ^> X r , or vice- versa. In the limit when 
A r ^> A s , the leading singular behaviour of the squared matrix element can be written as in 



Eq. (|6.2j), which coincides with Eq. ( [4.17 ). In the limit when A s ^> A r , the leading singular 



behaviour of the squared matrix element can be written as in Eq. ( |6.36| ) which is also the 
same as Eq. ( 4,17| ), 



There is yet another strongly-ordered soft-collinear limit that occurs when a soft gluon 
of momentum p^ = X s q makes a collinear splitting into a quark- antiquark or gluon pair 
of momenta pi and p r . The collinear splitting is parametrized as in Eq. ( |3.1[) with rescaling 
of k r ± as k r ± = X r k T ±. We consider the limit when A s — > with A s » A r . When the 
splitting of the gluon is into a qq pair, we start with the factorization formula for soft-qq 
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emission, Eq. (|4.20[ ), suitably re-labelled, 
sg*|-M^ 8 (Pi,JV,...)l !I = (6-60) 
= ^a a ^f\- £ E ( SjiSlr + SHSjr ~ SjlSir - 2 S -pA T R \M^ UJ0 (. ■ -)| 2 • 



j(ir) s l(ir) Sj(j r ) 



The strongly-ordered limit can be obtained by substituting into Eq. ( |6.60 ) the Sudakov 
parametrization of the momenta as in Eq. (|3.1|) , which leads to 



S( ir )C ir |A4i2+ 2 (pi,!*p, • • -)| 2 = 

= {8 7ra s ^) 2 5 M 5 frq j- £ ^MW^W)^^..^, (6.61) 

i.e. to the soft limit of the collinear factorization formula (Eq. ( |6.37| )). When the splitting of 
the gluon is into a gluon pair, we start with the factorization formula for doubly-soft gluon 



emission, Eq. ( 4.21 ), again suitably relabelled. The abelian part does not have a leading 
0[1/(A 2 A 2 )] singularity. Substituting the Sudakov parametrization of the momenta leads 



to an equation analogous to Eq. ( 6.61| ) with the pjj^ — > Pgg substitution. 



7. Singly-unresolved limits of the factorization formulae for doubly- 
unresolved emission 

We are now ready to list all the terms that are necessary for the explicit construction of 
the subtraction term A^l-M^i^l 2 - 

7.1 Singly-collinear limits of the doubly-unresolved factorization formulae 



First we consider the limit of Eq. (4^4) when partons i and r become collinear. In order 
to describe the ir — > i + r collinear splitting, we introduce the Sudakov parametrization 
of the parton momenta as in Eq. ( |6.45 ). The collinear region is identified by performing 



the rescaling k — > \n and studying the limit A — > 0. This limit is the strongly-ordered 
doubly-collinear limit, therefore, the factorization formula Cj rs |.M m +2| 2 factorizes as 

C ir C ira \M!® +2 (pi,Pr,Ps, • • -)| 2 = C( ir ) 8 Ci r \M!® +2 (pi,Pr,P 8 , ■ ■ -)| 2 , (7.1) 

where the right-hand side is the iterated collinear limit of the collinear factorization formula 
given in Eqs. (|6i3C|) - (|B3^) . 



We parametrize the collinear limit for the other doubly-unresolved factorization for- 



mulae as in Eq. (3.1). The factorization formula that is valid in the doubly-collinear limit 



(Eq. (|4.15|) ) is itself factorized in the two independent collinear splittings. Therefore, the 



limit of Eq. (|4.15 ) when parton i becomes collinear to r leaves that equation untouched, 
namely 

Ci r Ci r - js \M^ +2 {pi,Pr,Pj,p s , • • OP = GirCjs\M^ +2 (pi,p r ,Pj,p s , ■-.)?■ (7.2) 
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The same is true for Eq. ( 4,17] ), 



C ir CS ir . >s \M^ +2 (pi,Pr,Ps, ■ ■ 01* = CSir;s\M^> +2 (pi,Pr,Ps, ■ ■ ■ 



(0) 



(7.3) 



The collinear limit of the factorization formula for the emission of a soft gg-pair is 
clearly the strongly-ordered limit. The same is true for the non-abelian part of the two- 
gluon soft emission, thus, 



C^lM^p^Pr, ...)\ Z = S iir) C ir \M ( Z 2 ipi,Pr, ■ • 01* , 



(0) 



(7.4) 



where the right-hand side is the soft limit of the collinear-factorization formula given in 
Eq. ( |6.37 '). The abelian part of the two-gluon soft emission cancels in C/ ct A2|A / f^_ 2 | 2 > 
therefore, we do not record it. 



The computation of the remaining terms in Eq. ( 5.10| ) is straightforward. The action 
of C{ r on Eq. ( 4.32| ) reproduces the same, 



Cj r Cj r . s CSj r;s |.A4£ rl , 2 {pi,PriPsi • • 0| 2 



and so does it when applied on Eq. ( 4.36 ), 

CirCir-jgCSir-slM^J^iPij Pri Pj > Psj • • OP 

From Eq. fl4~37j ) we get 

C rs C irs S(,f\M^ +2 {pi,Pr,Ps,.. 



Ci rs CSi r - !S \A4^J +2 (pi, p r , p s , • • -)| 2 



C ir; j s CSi r - s \M^ +2 (pi,p r ,Pj,p s , . . .)\' z . (7.6) 



(7.5) 



SvrQ s/ u 2e ) 2 T 4 2 T R 



and from Eq. ( [4.41 ), with z r = 1 — z s , we find after some algebra 

C rs C irs St h) \M^ +2 ( Pi ,p r ,p s , . . 0| 2 = 
= (8vra s /i 2e ) 2 x 



\M^(p t ,...)\ 2 , (7.7) 



xT 2 C A 



2zi ( z r z s \ . . s ik 
— — + — + (1 - e) — 

Z r -(- Z s y Z s Z r J SrsSi(rs) 



\Mt ] {Pu ■■■)?■ (7.8) 



7.2 Singly-soft limits of the doubly-unresolved factorization formulae 

We consider the limit when the soft momentum is parametrized as p r = \q r and let A — > 
at fixed q r . In this limit the triple-parton splitting kernels factorize into a soft function 
times the Altarelli-Parisi kernel obtained from the triple-parton kernel by deleting the soft 
parton. Therefore, we obtain 



S r C i rs | M.U +2 (pi ,p r ,p 

,2e D (S) 



87TQ S /J, 2£ Pr i f (Zi, Z r , Z s , Si r , Sis, Srs'i Cj s | A4^J + ^ (pi,Ps, • • I* i 
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with 



P n n' n' (^*> ^Ti Sir i ^is i Srs ! e) — , 



2 % 



-^qifjrgs (^i j ? 5 ^ir 3 ^is 3 *^rs 3 ^ ) — ^F _ H~~ Ca 



2s, ; 



+ 



1 2= 1 Z,; 



Sfg Syg ™y S%y 



Pgig r g s %ri z si $iri $isi Sy S ] — C*A ( 



Sis 1 %i 1 z s . 

S °iySy§ S ' %y Zy Sy g Zy 



(7.10) 
(7.11) 



PgrlsQi ( Z ii z r, z si Sir, ^isi s rs'i e ) — Cf 1" C*A f 1 ) j (7-12) 

SiySy S \$rs z r S{y S{ySy S 



(7.13) 



The soft limit of the soft-collinear factorization formula in Eq. ( [4.17 ) leaves it un- 
changed if p s — > 0, 



S s CS ir . >s \M^ +2 (pi,p r ,Ps, ■■■)\ 2 = CSir;s\M^ +2 (pi,p r ,P 



(0) 



(7.14) 



However, it is crucial that in Eq. ( 7.14| ) we use Eq. ( |4.19| ). When only either parton r or s 
becomes soft, the asymptotic behaviour of the factorization formulae for the emission of a 
soft qq pair (in Eq. ( 4.2Q )) or of a gluon pair (in Eq. ( |4,21| )) is 



S r S r f\M { °UPr,Ps,---)\ 2 = 0, 

S r Si^\M^ +2 (p r ,p s , . . .)| 2 = S r S s \M^ +2 ( Pr , Ps , 



(7.15) 
(7.16) 



i.e., the iterated soft limit of Eq. (|6.13|). 



The computation of the remaining terms in Eq. fl5,9|) is straightforward. From Eq. ( 4.41 ) 
we get 

S s C irs S r T h) \M^ +2 (p uPr , Ps ,...)\ 2 = 
= (8vra s ^ 2e ) 2 x 



1 2z. 



1 2zi 



x Tf C A + 



1 2zf 



Sir Syc Zq Sj,o Syc Zy Sj_y Sic Zy Z, 



\M®(pi, ...)| 2 - (7-17) 



The action of S s on Eq. ( 4.34| ) leads to the same, 

(Pi,p r ,p s , ■ ■ •) 



and so does it when it is applied on Eq. (4.45), 



(pi,Pr,Ps, ■ ■ ■ )\ 2 = CS ir . iS Srs\M^ +2 ( Pi ,Pr,Ps,---)\ ■ 



(0) 



(7.18) 



(7.19) 



Finally, we consider the limit when the soft momentum r is also collinear to the mo- 
mentum of parton i. Clearly, in Eq. Q the operator C*. acts only on the pg> functions, 



Ci r S r Ci rs I A'l^^ (.Pii PriPsi 



2s 



)| C is \M^ +l {pi,p s ,...)Y . (7.20) 
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fs) 

In this limit, for the non- vanishing P- ■ ' functions we obtain 



Ci r P^ rs (^Zi, Z r , Z s , Si r , Si s , s rs , — ' . (7.21) 

Sir Z r 



The action of Cj s on Eq. ( 7,14| ) is 



Cj s S s CSi r - :S | A^^j^j_2 (Pi) Pj i Pr i Psi 

s^fc)3_L- ^ T 2 (^(ftr, . . 01 P$ p |A*g>(Pir, •••)>, (7-22) 



while on Eq. ( |7.16| ) it acts like 

C JS S r S r ^\M^ +2 (Pj,Pr,Ps,...)\ 2 = S r C JS S s \M^ +2 (p 3 , Pr ,p s ,...)\\ (7.23) 
which is given explicitly in Eq. ( 6.17| ). Furthermore, $ 



c rs s s s(r b) i^L(Pr,p s ,. 



2 



»^ 2 f^A^EE |^(r)|^ (j . (...)| 2 . (7.24) 



Applying the operator C JS S S on Eq. ( 4.45| ) we get 



Cj s S s CSj r;s S rs | A / I m ^2 i.Pit Pj i Pri Psi • • 01 



2e fT 2 J_ 2* 2 J_ 2^ |-M M (p . ^. _ _ >)|2 j (? _ 25) 



which holds also for j = i. Finally, the action of C rs on Eq. (fT-17 ) is 



C rs S s C irs S r T h) \M^ +2 ( Pl ,p r ,Ps, • • 0I 2 = {^a s ^fn C A — — \M$fa ^ 

Si r S rs z s 



(7.26) 



As we noted at the end of Sect. ||, the factorization formulae of this section, with similar 
ones that are not listed, can be used to write the terms on the right-hand side of Eq. (|5.£| ) 
in the form given in Eqs. (|5.9| )- (|5.11| ). 



8. The singly- and doubly-unresolved limits of the NNLO subtraction 
term 



We now turn to the proof of Eqs. (|5.2j)-( |577| ). The explicit check of these equations involves 
straightforward, but lengthy algebra. Therefore, we only outline the proof, which uses the 
fact that Ai2|A4^ 2 | 2 = Ai A2 1 A^^^_ 2 1 2 ? anc ^ ^ the interested readers go through the 
details. 



"'"In the C r sSs limit the momentum p r is the correct variable in the factorized matrix element (see 



Eq.(3.16)) 
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Let us first consider the singly-unresolved limits. Grouping the terms as in Eqs. fT2| ) 
and (|]§), we obtain 

X^A^M^l 2 = (x s A 2 + (XsA, - X S )A 2 ) \M%\ 2 , (8.1) 



where X s stands for either the collinear limit Cj 8 , or the soft one S s . As a result, 

X s (A 1 +A 2 -A 1 A 2 )\M^ l ) +2 \ 2 = Xsl-M^l 2 

+ (X s Ai - X s )(l - A 2 )|^2+ 2 | 2 • (8-2) 

The factor in the first parenthesis on the right-hand side contains operators that take the 
iterated singly-unresolved limits, which give vanishing contributions on (1 — A 2 )|.M^ 2 | 2 , 

(XAx - X s )(l - A 2 )|_M^ 2 | 2 = • (8.3) 

We checked Eq. flS.3| ) explicitly for all the terms in Eqs. ( |3.21 ) and ( |3.22| ). It is a conse- 
quence of Eqs. (|4.59|) -( 4.62 ), which hold also in the strongly-ordered doubly-unresolved (or 
equivalently in the iterated singly-unresolved) regions. Thus only the first term on the 
right-hand side of Eq. ( |8.2|) remains, and that is the relevant singly-unresolved limit of the 
squared matrix element. Therefore, Eqs. ( |5.2[) and (|5.3|) hold. 



Turning to the doubly-unresolved limits, we first observe that according to Eqs. ( |4.59 )- 



(|4.62[) we have 



Y rs A 2 \M { Z 2 \ 2 = Y rs \M { ° ) +2 \ 2 > (8-4) 



where Y rs stands for any of the operators Cj rs , Ci r -j s , CSi r - s and S rs . Thus, all we have 
to prove is the following relation: 

Y re (A 1 -A i3 )\M®. 2 \ 2 = 0. (8.5) 



We checked the validity of Eq. ( |8.5| ) explicitly for the four different limits. In general it 
follows from 

F rs (Ai - A^IA^I 2 = Y rs A x (l - A 2 )|A4^ 2 | 2 , (8.6) 

and observing that in any term X s |.A/f^+ 2 | 2 in the NLO subtraction Ai|.M^ +2 | 2 the de- 
pendence on s is already factorized, therefore, the terms Y rs X s\M^ +2 \ 2 are factorization 



formulae valid in the iterated singly-unresolved limit, or according to the results of Sect. |6.5| , 
in the strongly-ordered doubly-unresolved regions. In these regions Eqs. ( 4.59| )-( |4,62; ) hold, 
therefore, 

F rs A 1 (l-A 2 )|>f^ 2 | 2 = 0, (8.7) 



which leads to Eq. 

It is worth noting that in Eqs. (|8.3| ) and fl3.5|) we have to take unresolved limits of 
factorization formulae that are valid in different unresolved regions of the phase space. For 
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instance, Eq. ( ]8.3| ) contains terms such as S s C sr S r A 2 |>I^ 2 | 2 . According to Eq. (gup 
such terms contain the factorization formula C S rS r Ci rs \A4^ +2 \ 2 1 given explicitly by 

C sr S r .Cf rs |A / (^ n |(_2(pi,Pr)Ps) • • -)l = 

= (8vra s /i 2e ) 2 x 

<Af#(ft., • • -)\P^M, k s± ;e)\M^(p is , ...)}. (8.8) 

S rs Z r Si s 

In the region where p s — > 0, this expression contains the ambiguous limit lim Zs ^o z s/ z s- 
This limit can be computed using the constraints the momentum fractions have to obey. 
The z s in the numerator is relevant in the collinear limit of partons s and r, where it fulfills 
the constraint z s + z r = 1. The z s in the denominator is relevant in the region where the 
partons i and s are collinear when z s + %i = 1. Thus the momentum fractions, needed to 
define the various collinear limits, have to be extended to other parts of the phase space 
such that they fulfill different constraints in different regions. For instance, using the energy 
fractions, mentioned at the end of Sect. 4.1, we find that the soft limit of Eq. ( |S.§| ) is 

S s C sr S r C irs \M^ +2 ( Pi ,Pr,Ps, ...)\ 2 = (8vroV £ ) 2 ^ C K Tl\M^{pi, . . .)| 2 . (8.9) 



S r aSi 



The cancellations discussed in Eqs. Q8.3] ) and (8.5) all take place with momentum fractions 
defined as energy fractions. 



9. A possible general subtraction method for computing NNLO correc- 
tions to jet cross sections 

Building on our studies in the previous sections concerning the matching of the various 
singly- and doubly- unresolved limits of the squared matrix element, in this section we 
outline a possible general subtraction method for computing NNLO corrections to jet cross 
sections. 

The full cross section to NNLO accuracy can schematically be written as 

^LO i _NLO , NNLO 1 \ 

<7NNLO = CT +<J +a . (9.1) 

Assuming an m-jet quantity, the leading-order contribution is the integral of the fully 
differential Born cross section dcr^ of m final-state partons over the available m-parton 
phase space defined by the jet function J m , 

<J LO = [ d<j*J m . (9.2) 

J m 

The NLO contribution is a sum of two terms, the real and virtual corrections, 

^NLO _ / dcr^ +1 J m +i + / da^Jm ■ (9.3) 

Jm+l Jm 
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Here the notation for the integrals indicate that the real correction involves m + 1 final- 
state partons, one of those being unresolved, while the virtual correction has m-parton 
kinematics. The NNLO correction is a sum of three contributions, the double-real, the 
one-loop single-unresolved real-virtual and the two-loop double-virtual terms, 

J m+ i + / da^J m . (9.4) 

Jm+2 J m+1 J m 

Here the notation for the integrals indicate that the double-real corrections involve m + 2 
final-state partons, the real-virtual contribution involves m + 1 final-state partons and the 
double virtual term is an integral over the phase space of m partons, and the phase spaces 
are restricted by the corresponding jet functions J n that define the physical quantity. 

The various contributions to the NLO and NNLO corrections are divergent. We regu- 
larize these singularities using dimensional regularization, that is defining the integrals in 
d = 4 — 2e dimiensions. After ultraviolet (UV) renormalization of the one- and two-loop 
matrix elements involved in dcr v , dcr RV and dcr vv all the contributions to cnnlo are UV 
finite. In general, the jet function of m + n partons vanishes if n + 1 or more partons 
become unresolved. Explicitly, J m vanishes when one parton becomes soft or collinear to 
another one, 

J m (pi, . . . ,p m ) -»• , if pi-pj^O, (9.5) 
and Jm+i vanishes when two partons become simultaneously soft and/or collinear, 

J m+1 (px, ■ ■ ■ ,p m +i) -> , if Pi -Pj and p k -pi -> (i ^ k) . (9.6) 

Therefore, the Born contribution du^ is integrable over the one- or two-parton infrared 
(IR) regions of the phase space. Also the real and virtual contributions d<7^ +1 J m+ i and 
da^Jm vanish when we integrate over the two-parton IR region of the phase space, but are 
separately divergent in d = 4 dimensions when only one (real or virtual) parton becomes soft 
or collinear to another parton. Their sum is finite for infrared-safe observables. Formally, 
infrared safety is expressed by the following requirements for the jet function, 

Jn+l(pi, ■ ■ ■ ,Pj = M, ■ ■ ■ ,Pn+l) -»• JniPl, ■ ■ ■ ,Pn+l) if A^O, (9.7) 

J n +l(pi, ■ ■ ■ ,Pi, ...,Pj,...,p n +l) JniPl, ■ ■ ■ ,P, ■ ■ ■ ,Pn+l) if Pi -> zp , Pj ~> (1 ~ z)P 

(9.8) 

for all n > m. Equations ( |9.7D and Q9.8|) respectively guarantee that the jet observable 
is infrared and collinear safe for any number n of final-state partons, i.e. to any order in 
QCD perturbation theory if there is only one infrared parton. The n-parton jet function J n 
on the right-hand side of Eq. ( |9.7p is obtained from the original J n +i by removing the soft 
parton pj, and that on the right-hand side of Eq. ( |9.8[ ) by replacing the collinear partons 
{PiiPj} by p = Pi +Pj- For infrared-safe observables, process- and observable-independent 
techniques exist to cancel the IR divergences by devising an approximate fully differential 
cross section over the m-parton phase space, der^' l such that (i) dcr^'j^ matches the 
pointwise singular behaviour of dcr R , 1 in the one-parton IR regions of the phase space in 
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any dimensions (ii) and it can be integrated analytically over the one-parton phase space, 
so we can combine it with d<r^ before integration. We then write 



NLO 



m+l 



A R ' A T 
aa m+l J m 



+ 



R,A 



J n 



(9.9) 



where both integrals on the right-hand side are finite in d = 4 dimensions (in order to prove 
the finiteness of the first integral, we have to use Eqs. (|9.7| ) and ( |9.8[) ). The final result is 
that we were able to rewrite the two NLO contributions in Eq. fl9.3|) as a sum of two finite 
integrals, 

- NL ° " f doffiS + I da^° , (9.10) 



a 



m+l 



that are integrable in four dimensions using standard numerical techniques. 

The three contributions to the NNLO corrections in Eq. ( p.4j ) are separately divergent 
in both the one- and the two-parton IR regions, but their sum is finite for infrared-safe 
observables. To NNLO accuracy, the formal requirement of infrared safety is expressed by 
the following relations, 



-2(P1, ■ 
-2(P1, ■ 
if 

-2(Pl, ■ 
if 



• ,Pj = M, ■ ■ ■ ,Pk = Ar, . . . ,p n+2 ) -> Jni.Pi, • ■ 

• ,Pj = Xq,-.. ,Pk, ...,Pl,... ,Pn+2) ~> Jnipi, 

A -> and p k -> z k p k i , Pi ->• (1 - Zfc)PfcJ , 



Jn+2(Pl 



■ ■ iPii ■ ■ ■ iPj i ■ ■ ■ iPki ■ 
Pi -> ^iPij , Pj -> (1 
• • iPii ■ ■ ■ iPji ■ ■ ■ iPki ■ 



,Pl, ■ ■ ■ ,Pn+2) JniPU ■ ■ 

Zi)pij and p k -> z k p k i 

,Pn+2) — * Jn{Pl, ■ ■ ■ ,Pijk, 



:Pn+2) 
■ ,Pkh ■ 

■ ■ ■ iPiji 
Pi ~i 



if A — 

• ,Pn+2) 



0, (9.11) 



(9.12) 

• • ,Pkl, ■ ■ ■ ,Pn+2) 

[1 - z k )pki , (9.13) 



if pi -> ZiPijk , Pj zjpijk and p k (1 



• • ,Pn+2) 
Zj )Pijk 



(9.14) 



for all n > m (as implicitly mentioned before, J m +2 vanishes if three or more partons 
become unresolved). Eqs. ( 9.11 )-( |9.14| ) guarantee that the jet observable is infrared safe 
for any number n of final-state partons, i.e. to any order in QCD perturbation theory if 
there are exactly two unresolved partons. The n-parton jet function J n on the right-hand 
side of Eq. ( |9.11 ) is obtained from the original J n +2 by removing the soft partons pj and 
Pk, that on the right-hand side of Eq. ( |9.12| ) by removing the soft parton pj and replacing 
the collinear partons {p k ,pi} by pki = p k +Pu that on the right-hand side of Eq. ( 9.13; ) by 
replacing the collinear partons {pi,Pj} by pij = pi +pj, and {pk,Pi} by p M = Pk + Pi and 
that on the right-hand side of Eq. ( |9.14 ) by replacing the collinear partons {pi,Pj,Pk} by 
Pijk =Pi +Pj +Pk- 

In order to cancel the one- and two-parton singularities between the double-real, real- 
virtual and double- virtual terms in Eq. (p.4|) , we proceed like in Eq. (|9.9| ) , namely, we 



introduce approximate cross sections dcr^ + ' 2 2 and do",^'^ 1 and rewrite Eq. ( |9.4[) as 



NNLO 



+ 



m+2 

da vv 



aa J m +2 



i R,Ft,A2 t 
aa m+2 J rr 



+ 



da Jm+1 



m+2 



j RR,A2 j | 
aa m+2 J m + 



m+l 



m+l 

RV 
m+l 



, RV.Ai j 
aa m+l J n 



+ 



j RV.Aj T 
_i_i Jm 



(9.15) 
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The fully differential approximate cross section to the double-real emission is constructed 
such that it has the same pointwise singular behaviour in d dimensions as da^^ 2 i n the two- 
parton infrared regions (hence the subscript 2 on A in the superscript). Then Eqs. ( |9. 11 ) — 



(9.14) guarantee that the (m + 2)-parton integral in Eq. ( |9,15| ) is finite over the doubly- 



unresolved regions of the phase space. However, this integral is still divergent in the singly- 
unresolved regions, because both the fully differential cross sections d<7 RR and dcr^'^ 2 as 
well as the jet functions J m +2 and J m have different limits in the singly-unresolved regions. 

The real-virtual contribution has two types of singularities: (i) explicit e poles in the 
loop amplitude and (ii) singular behaviour in the singly-unresolved regions of the phase 
space. To regularize the latter, the fully differential approximate cross section, du^f 1 , 
to the real-virtual contribution must have the same pointwise singular behaviour in d 
dimensions as d<T R l ^ 1 itself in the one-parton infrared regions (hence the subscript 1 on A 



is 



in the superscript). Using such a regulator, the difference 
still not integrable over the full (m + l)-parton phase space in d = 4 dimensions, simply 
because the integrand still contains e poles in phase space regions away from the one-parton 
unresolved infrared regions. The pole structures of da™ +1 and dcr^X^ 1 are different. 

In order to cancel the remaining singularities in the (m + 2)-parton and {m + 1)- 
parton integrals in the first line of Eq. ( |9.15| ), we perform a subtraction procedure between 
these two integrals, which is very similar to that used in a NLO computation, Eq. (|9.SQ. 
This amounts to devising fully differential approximate cross sections to both dcr RR _ 2 and 
d fJ m+2^ 2 that match the singular behaviour of these cross sections in the singly-unresolved 
infrared regions in d dimensions, which are denoted by du^'^ 1 and dcr^'^ 12 , respectively. 
Thus, the full NNLO cross section is written as 



^.NNLO 



m+2 



^i RR t ^i RR,A2 j i RR,A} j , i RR,A-^2 j 

dC Jm+2 — d Cm+2 J m ~ a<J m+2 J m+1 + aa m +2 J m 



+ 



I m+2 L 



771+2 ,J Tfl 



+ 



+ / da RV J m+1 - da^'f 1 J m + / 

J m+l J m 

+ f da VY J m + f d^' 2 A V m + f da^KJ m . (9.16) 

Jm Jm+2 Jm+1 

The first integral on the right-hand side of this equation is finite in d = 4 dimensions by 
construction. Hence, this first integral can be performed numerically in four dimensions 
using standard Monte Carlo techniques. 

The singularities are associated to the integrals in the second and third lines of 
Eq. ( |9.16D . We shall show in a separate publication that the integration of [d<r m+ ' 2 1 J m+ i — 
da^^2 12 Jm] can be performed analytically over the one-parton phase space of the unre- 
solved parton of the (m + 2)-parton phase space, leading to e poles, which cancel exactly 
the remaining e poles in the integrand [da^_ 1 J rn+ i — da^^f 1 J m ]. As a result the integrals 
in the second line of Eq. ( |9.16| ) can be combined into a single (m + l)-parton integral, where 
the integrand is free of e poles and can be integrated numerically in four dimensions. We 
emphasize that this statement does not follow from unitarity, but has to be proven. 
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Once we have shown the finiteness of the combination of integrals in the second line of 
Eq. ( |9,16 ), it follows from the Kinoshita-Lee-Nauenberg theorem that the combination of 
the integrals in the third line must be finite, though separately all three are divergent. If one 
is able to carry out analytically the integration of da*^^ 2 over the two-parton subspace 
and that of dcr^'^ 1 over the one-parton subspace of the doubly- and singly-unresolved 
regions of the (m + 2)-parton and (m + l)-parton phase spaces, respectively, leading to 
e poles, we can combine these poles with those in the double-virtual contribution. As a 
result all the divergences cancel and the remaining m-parton integral can be performed 
numerically in four dimensions. 



The final result of the above manipulations is that we rewrite Eq. ( |9.4[ ) as 



a 



NNLO 



m+2 



i NNLO , 
do m+2 + 



da. 



m+l 



NNLO 
m+l 



+ / da 



NNLO 



that is, a sum of three integrals, 



da, 



NNLO 
m+2 



i RR t ^i RR,A 2 j i RR,A} j , n RR,A-^2 j 

U-°'m+2 J m+2 ~ ^ m +2 J m ~ m +2 J m+1 + aa m+2 J n 



£=0 



i NNLO 
QCJ m+l 



RR.Ai 



j RV j i RV,A^ j j I i RR,A} j j i\.i\,j-i^2 j 

fl(7 m+l J m+l — ^m+l + / [ CtcJ m+2 ^m+1 — clcr m+2 



£=0 



and 



da, 



NNLO 



, \ \' f n RR,A 2 . f , RV.Ai 
I" J 2 dc7 m+2 + _/ d<7 m+l 



<^m ; 



(9.17) 

(9.18) 
(9.19) 

(9.20) 



e=0 



each integrable in four dimensions using standard numerical techniques. 



The factorization formulae presented in the previous sections can be used in con- 
structing the approximate cross sections do"^'^ 1 , dcr^'^ 2 and dcr^'^ 12 . However, we 
emphasize that the momenta in the matrix elements on the right hand sides of the factor- 
ization formulae are unambigously defined only in the strict unresolved limits. Therefore 
the factorization formulae cannot be directly used as true subtraction formulae, unless 
we can explicitly specify which partons become unresolved — either by partitioning the 
squared matrix element as in Ref . , or using the measurement function as in Ref . Q , or 
by partitioning the phase space as in Ref. Q . These procedures lead to the so-called residue 
methods which are well understood in NLO computations. If one wishes to avoid specifying 
the unresolved partons, then one has to implement exact factorization of the phase space, 
such that the integrations over the unresolved momenta can be carried out and momentum 
conservation is maintained as done in the case of dipole subtraction method of Ref. §, 
which is worked out completely for NLO computations. 



10. Conclusions 



In this paper we studied the infrared structure of the known factorization formulae for tree- 
level QCD squared matrix elements in all the possible soft and collinear limits. We pre- 
sented new factorization formulae for the colour-correlated and the spin-correlated squared 



matrix elements that we termed iterated singly-unresolved limits. We pointed out that 
soft factorization formulae do not exist for the simultaneously spin- and colour-correlated 
squared matrix elements which indicates that in the general subtraction scheme envisaged 
by us, the azimuthally correlated singly-collinear subtraction terms must not contain colour 
correlations. This can be achieved naturally for those processes when the colour charges in 
the colour-correlated squared matrix elements can be factorized completely (see Appendix 
A). We derived the factorization formulae of the squared matrix element that are valid in 
the strongly-ordered doubly-unresolved regions of the phase space and showed that these 
are equal to the iterated singly-unresolved factorization formulae. We dealt only with fi- 
nal state singularities. Using crossing symmetry, it is possible to extend our results to 
cases when initial state partons become unresolved, which we shall present in a separate 
publication. 

We introduced a new formal notation for denoting the various limits in the overlapping 
unresolved parts of the phase space. Using this notation, we explicitly constructed subtrac- 
tion terms that regularize the kinematical singularities of the squared matrix element in 
all singly- and doubly-unresolved parts of the phase space and demonstrated that the sub- 
traction terms avoid all possible double and triple subtractions. The relevant subtraction 
terms were presented in Sects. ||| and [|. As a result the regularized squared matrix element 
is integrable over all phase space regions where at most two partons become unresolved. 
We emphasize however, that the subtraction terms derived in this paper can be used as 
true subtraction terms only if the unresolved momenta are specified, like in Refs. |4| in 
the case of NLO computations. In an approach, where any parton can become unresolved, 
the factorization formulae have to be extended over the whole phase space. This extension 
requires a phase-space factorization that maintains momentum conservation exactly, as 
done in Ref. || in the case of NLO computations, but such that in addition it respects the 
structure of the delicate cancellations among the various subtraction terms. So far none of 
these methods has been demonstrated to be applicable in practical computations of NNLO 
corrections. The kinematics of two-jet production in electron-positron annihilation is suffi- 
ciently simple so that we were able to implement the subtraction scheme presented in this 
paper for the case of m = 2. We computed the cross sections in Eqs. ( |9.18| ) and ( 9.1S| ) and 



found numerically stable results for any observable. The details of the computations will 
be given elsewhere. 
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A. Explicit computation of the soft limits of known colour-correlated 
squared matrix elements 



In this appendix we check the validity of the soft factorization formulae Eqs. (p.cf ) and (|6.1l| ) 
using the colour-correlated squared matrix elements for the e + e~ — > qq + ng processes 
(n = l, 2). 



A.l Soft limit of \M 



(0) | 
3;(i,fc)l 



Using colour conservation, we can compute the product of the colour charges acting on the 
three-parton colour space in terms of Casimir invariants ||, 

J-q-J-g — ~2 °F' 1 9 ' 1 q — 1 9 ' 1 q ~ ■ 

Thus the colour charges in the colour-correlated squared matrix elements for the e + e~ 



(A.l) 



qqg process can be factorized completely [38]. If the quark, antiquark and gluon is labelled 
(in this order) 1, 2 and 3, then we have 



(0)|2 
3 I 



and 



I 2 _ \ M (0) ,2 
3)1 — l- /vt 3;(2,3)l 



— lM 0) ! 2 



(A.2) 



(A.3) 



Consequently, the limit of the colour-correlated squared matrix elements when p% gets soft 
is simply colour factors times the soft limit of the squared matrix element [39], 

I^ 1>2) (P1,P2,P3)| 2 - -8na s ^ -Cf) S 12 (3)\M { ^ 2) (pi,p 2 )| 2 , (A.4) 

which — using T\ ■ T2\M 2 °\pi,P2)) = — Cf\ M. ® {PI1P2)) — can also be obtained from 
Eq. (|6ij|), and 



M^ipi^Ps)? - 87ra s /^ ^S l2 (3)\M^ 2) (pi,p 2 )\ 



C A 



,(o) 



(A.5) 



which can also be obtained from Eq. (6.11). 



A.2 Soft limit of \M^ ] (ik) \ 2 



The colour-correlated squared matrix elements for the e + e~ — > qi&gzdi process were com- 
puted in Appendix B of Ref. pQV We rewrite the expressions presented there using our 
notation: 



M 



(0) 

4;(i,fc) 



1 



N C C F ( C| M* fe + C A C F M l x k + C\M% ) . 



(A.6) 



§ Note a typographic error in Eqns. (B11-B13) of Ref. J4(J: the labelling convention given in that paper is 
correct for the process e + e~ — + qxqigzQA we consider here. In order to get the Mg k , M lk and MJ* matrices 
relevant to the labelling convention used in Ref. ^(J the 2, 3 and 4 indices of the M^ k matrices should be 
cyclicly permuted, (2,3,4) -> (4,2,3). 
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where the non- vanishing elements of the matrices M$ k , M x k , M xx are given by 



M 12 = 2|m 3 | 2 , 



(A.7) 



M 



12 



3\m s \' 2 + \ mi \ 2 + \m 2 \' z , M x 13 = M x i4 = Mf = Mf = |m 3 | a , (A. 



r 23 



r 24 



34 



-(|mi| 2 + |m 2 | 2 ) 



M x 12 = |m 3 | 2 -M^ 



(A.9) 



Mil = Mil = -\(\m?,\ 2 - Irml 2 ) , M l x l = Mf x = ~(\m 3 \ 2 - |m 2 | 2 ) , (A.10) 

and the |mj| 2 functions are the helicity-summed squared helicity amplitudes with all cou- 
pling factors included. In terms of the amplitudes defined in Appendix A of Ref. [ 40 1 

|mi| 2 = ^|m(l,3,4,2)| 2 , \m 2 \ 2 = |m(l, 4, 3, 2)| 2 , (A.ll) 



hel 



hel 



and m 3 = mi + m 2 - Using the soft factorization formula for the helicity amplitudes, it 
is not difficult to obtain the following factorization formulae for the \rrii\ 2 functions when 
p 4 -»0: 



\ r mi(pi,p 2 ,P3,P4:)\ 2 ^ 8vra s ^ 2E 5 23 (4)— — \Mf ] {p 1 ,p 2 ,p-i)\ 2 , 

1\ c Cf 

\m 2 (pi,P2,P3,P4)\ 2 k 8vra s ^ 2E 5i 3 (4)^-^-|A^^ 0) (pi,p 2 ,P3)| 2 , 



\m 3 (pi,p 2 ,p3,p4)\ 2 ~ 87ra s ^ 2£ 5i2(4)-^^-|A4^ 0) (pi,p 2 ,P3)! 2 • 



(A.12) 
(A.13) 
(A.14) 



Substituting Eqs. $7?\)- QKI$ ) into Eq. < KM , we obtain the soft factorization formulae for 
the colour-correlated squared matrix elements, when p 4 — > 0: 



^4;(i,ifc)(Pl>:P2,P3,.P4) 



-8-7ra s /i 



2£ 



■M3 (P1,P2,P3) Fik(pi,P2,P3,P4) , (A.15) 



where 



•^12(P1,P2 


P3 


Pa) 


= (c F 


^13(Pl,P2 


P3 


Pa) 


_ C A " 
2 


^14(P1,P2 


P3 


Pa) 


_ Ca ' 


J 7 34(P1,P2 


P3 


Pa) 


-h 



Ca 
2 



(C F - C A )S l2 (4) + ^ («S 13 (4) + 623(4)) 



5 12 (4) + ^5 23 (4) 
C F -^) S 12 (4) + ^S 13 (4) 



(A.16) 
(A.17) 
(A.18) 
(A.19) 



Using Eq. ([All) it is straightforward to check that the results in Eq. ( |A.15[) , Eqs. ( |A. 1€| ) 
and (lAl7|) can also be obtained from Eq. (|6.8| ), and those in Eqs. ( |A.1E| ), ( |A.18[) and 
( A.lSj ) can also be obtained from Eq. ( |6.11| ). 
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B. Collinear limit of the spin-polarisation tensor for the process e + e 



qqg 



Consider the tree-level matrix element for the production of a quark-antiquark pair and a 
gluon in electron-positron annihilation, 



Ml\: c i^ sl ' S2 ^(pi,P2,p) = (cic 2 c 3 | ® (s 1 s^ 3 \\Mf\p 1 ,p 2 ,p)} , 



(B.l) 



where {01,02,03}, {si, 82,^3} and {51,92,53} are the colour, spin and flavour indices of the 
partons. We define the spin-polarization tensor as in Eq. (|6.20|): 



%lsg(PUP2,P) = {M < 3 0) {p 1 ,p 2 ,p)\p){v\M 3 U> {pi,P2,p)) 



(0), 



(B.2) 



so -g^T q %( Pl , P 2,p) = \Mf\ Pl , P 2,p)\ 2 . 

In this Appendix, we examine the collinear limit of (fj,\Pgg\z, k±; £)\v) %fq^(pi,P2,p) 
(we suppress the second argument of the Altarelli-Parisi kernel, which equals 1 — z). From 
Appendix D of Ref. || we have (in d = 4 dimensions) 



%^o(Pl,P2,P) 



1 



yfp + vl v + 2 yi2 

where yij = Sij/Q 2 = Sy/(pi + p 2 + p) 2 and 



Mf\ Pl ,p2,p)\ 2 T^, 



(B.3) 



rpflV _ 1 ^2 _|_ 2^2^1 



+ 



g 2 Q 2 yi P Q 2 
2/12 - (yi2 + y2 P ) 2 



2 mpPiPi_ _ 2 mpP2P2 + 



yi P 



PiP u p^Pi 



V2 P Q 2 

yi2 - (yi2 + yip) 2 



Q 2 



Q 2 



+ 



V2p 



P2P U , £^2 

g 2 g 2 



+ 



1 



+7j (y?p + y^) 5' 



(B.4) 



We also have (in d = 4 dimensions) 



I « >(0)x x,2 n 8lTas + y 2p + 2 2/l2 I . ,(0),2 

I-m 3 (pi,P2,p)r = cv — ^ — - — \Mh ;| - 



g 2 



yipV2p 



Here l-M^I 2 is the squared matrix element for the production of a quark-antiquark pair 
in electron-positron annihilation, averaged over event orientation (so it has no dependnece 
on parton momenta). Substituting Eq. ( B.5 ) into Eq. (|B.3| ), we find: 



(B.5) 



T q %(pi,P2,p) = -8vra s C F -L \Mf\ 2 

yipy2p 



rj-i(J,U 



(B.6) 



The tree-level Altarelli-Parisi splitting kernel for g — > gg splitting, (/j,\Pgg(£, k; e)\v) 
is given in Eq. ( |3.5[) . The case of interest to us is when p in the spin-polarisation tensor is 
the collinear direction of the 5 — > gg splitting. In particular this means that p ■ k = 0. 
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With the explicit expressions for T^ u g (pi,p 2 ,p) and {fi\Pgg ((, k; £)\v), we can carry out 
the contraction 



2C A {\Mf\ Pl ,p 2 ,p)f 



i-C c 



(B.7) 



+87ra s C F \M 2 



(0)|2 



1 -£ 

S34 



SlK s 2k 



S34 



— - — ) +7iC(l-C) 
sip s 2p J Q 2 



SlpS 2p 



where we have introduced the following notation: 

Sin = 2p{K, S{p = 2piP, I = 1, 2 . 



S34 



C(i-C) 



(Bi 



We now derive the pi | |j? collinear limit of Eq. flB.?] ) , which is precisely defined by the 
following Sudakov parametrization of the momenta: 



£.2 h.2 
^ = Zl pM + ^-^jL_, = (1 _ ^) P M _ fcM 



Sip 



0l(l - Zi) 



21 2Pn 
= 0(4), fc ± ^0, 



1 - z\ 2Pn ' 



(B.9) 



where Pk± = nk± = P 2 = n 2 = 0. When calculating s\ K , we use pK = to eliminate the 
Pk term from s\ K . Then we get 



SlK 



-k± ■ k + 2 



21 



1 \ TIK 



Then the collinear limit of Eq. ( |B.7| ) reads 
C V p{f,\P$((,K-,e)\v)T q f g (p u p 2 ,p) = 
= 8ira s \M { 2 0) (P,p 2 )\ 2 2C A —x 

Sip 



k\ = 0(k A 



(B.10) 



(B.ll) 



1 - C C / «34 



-^ + S3 4C(i-C)(i-2i) 

Sip 



We can compute the same limit using the expression derived for the collinear limit of 



the spin-polarisation tensor in Sect. 3J]. Defining the pi\\p limit as in Eq. ( |B.9| ), we have 
(see Eq. (1^22|) ): 



T q f g (pi,P2,p) * 8na sf i 2£ — (M { 2 0) (P,p 2 )\ P% (z u k ±] P,n,e) \M 2 "> (P,p 2 )), (B.12) 

Sip 



where (r\Pgg((,K;p,n,e)\s) is given in Eq. ( 6.3l| ). Using Eq. ( [B.12] ) we obtain for the 
collinear limit of Eq. (IT~ 



Cip (»\P$((,K;e)\v) T^( Pl ,p 2 ,p) 



qqg^ 
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87ra s/U 2£ |^ 0) (P,p 2 )| 2 2C A — x 

Sip 



(B.13) 



l-C c 

1-e 



+ 2C F 



S34 



l-si 2 ^zi{k±^Kf (Pn)(nn) 

-K —2- — \-(l — Zi) 



1-zi 



Pn 



To see that Eqs. ( |B.11| ) and ( |B.13| ) are in fact equal up to subleading terms, note the 
following. From pn = 0, we get 

Zlk " ™ — * (B.14) 



Pk 



l-zi 



k± ■ k + 



(l- Zl ) 2 2Pn 



0(k A 



so the last term in the square brackets in Eq. ( B.13j ) is subleading. Furthermore, we see 
from Eqs. and flB.10|) that 



(fc± • «) s 
k 2 . 



1 ~ S L 
4zi sip 



+ 0(fcj 



(B.15) 



Using Eqs. ( gig ) and ( p7T5|) in Eq. ( pl3|) we obtain Eq. pH^ . This gives a direct check 
of result Eq. (6.31) for the splitting tensor. 



C. The collinear limit of the spin-polarization tensors using a helicity basis 

In this appendix we record the collinear limit of the T hh ' tensors using a helicity basis. 
There are two ways to find the factorization formula. One is to work directly with helicity 
amplitudes and use the known factorization formulae with the splitting functions factorized. 



The other is to project the factorization formulae of Sect. 6.3 onto a helicity basis. We 
explore both ways and show that these lead to the same expressions. The splitting functions 
are defined in d = 4 dimensions, therefore, in this appendix we set e = everywhere. 

Let us start with the former. In our notation the helicity amplitudes are obtained as 

4T(1 A \ ■ • • , ™ Am ) = (Ai, • • • , X m \Mt\pu ■ ■ ■ ,Pm)) ■ (C.I) 

The helicity amplitudes can be written as a sum of products of colour factors C a and 
partial amplitudes ^£+1, 

AT (l Al , • • • , ™ Am ) = 9T~ 2 E C a({l,...,m})< ee (^« lAl > • • • > ») > ( C -2) 

where <t({1, . . . , m}) denotes a certain permutation of the labels. The allowed permutations 
are not important in our present calculation. 

The collinear-factorization formula for the partial amplitudes when partons i and r 
become collinear reads (see e.g. Rcf. JO]]) 



^Si(iV*V • •) = Split te A e (^, r x r)AT(pi, •••) + ■••, 



(C3) 
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where the ellipses mean subleading terms, which are neglected in the collinear limit. Sub- 
stituting Eq. ( p.3| ) into the definition of the helicity-dependent tensor 

(C.4) 



T%&{pi,Pr, ■■■) = (M^ +1 ( Pi , Pr , . . (\' i \M < $ +1 (pi,p r , • • •)) 



we obtain the factorization formula 



C ir r^l( Pi ,p r , ...) = 8na s T 2 r £<A|P^(* r , (ir))\\')AT(pL ■ • •) ^T(p^, ■ ■ ■] 



xx 1 



where 



(\\P;£(z r ,(ir))\\') = - Spht^^S A) Spht^^sr^; 



(C.5) 



(C.6) 



is the representation of the Pj^ splitting tensors on a helicity basis and (ir) is the spinor 

product as defined, e.g., in Ref. [El]]. Using the known splitting functions for the processes 

~ A A' * A • A' 

q — * q r +gi and g — > g r +g r , we can easily compute the tensors (A|P 9i q r l |A ) and (A|P 9 / gr l |A'). 
For the case of quark splitting we have 



(A|^|A') = 5aa^#, 
therefore, we have to compute only the following cases, 

1 1+z 2 



(C.7) 



P"M-= p— 



P+~ = I P 

9i1r \ 9i1r 



1 ( Z l-Z 
+ — =— 



* 1 z 



For the gluon splitting we find 



(+\P++(z r ,(ir)) 
(+\P+-(z r ,(ir)) 
(+\P g --(z r ,(ir)) 

(+\p- t t(z r ,(ir)) 
(+\P+; r (z r ,(ir)) 



+ 



+ 



+ 



( + \P++( Zr ,(ir})\-y. 



+ 



(ir) 2 1 — z 

\Pg~g r (Zr,(i r )) 

\Petg~MAir)) 
\P+-{z r ^r)) 

\ P g~g r ( z r, (ir)) 
\Pg~g r (Zr,{ir)) 



(C.8) 



+ 



+ 



+ 



1 z 6 



1 ' 
Sir -L % 

o, 

1 f(l-z) 3 



+ 



1 



z z(l — z) 



2 l-z 



Sir % 

( + \P~l( Zr ,(ir))\+y 
1 z 
(ir) 2 1 — z ' 

<-|^K(zr,<<r»|+>- 
z(l-z). 



(ir)- 



(C.9) 



We can derive the factorization formula in Eq. ( |C.5| ) from Eq. ( |6.22| ) by inserting the 
resolution of the identity in terms of projection operators, 



5>><ai = i, 



(CIO) 
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and projecting the Lorentz indices onto helicity states with polarization vectors = (A|/i), 



Ci r T m+ l(pi,p r , • • •) 



8tt« s — J2(M$(p ir ,...)\\) (A|//) x 

Sir AA' 

x(A l |«) ( M |^5 p (« r ,Avj.;Pir,nir)>H </?|A$> HA') (\'\M$(Pir, ■ ■ •)) • (Cll) 



Comparing it to Eq. ( C.5 ), we see that 



Tl (X\P^}(z r , (ir))\\>) = — (A| M )(Ai|a) {^ fr {^ k r ±; P ir,n ir ))\u) (f3\\>) HA') . 

Sir 

(C.12) 

In order to compute the right-hand side of the equation above from Eqs. (|6.31| ) and (|6.32| ), 
we note first that the gauge terms do not contribute. Secondly, in Eq. ( |6.32| ) we left out 
terms that do not contribute to either Eqs. ( |6.25 ) and ( |6.26 ) or Eq. ( 6.33; ). Including these 
terms and neglecting the gauge ones, we have 



(s|Af (z,kr,P,n ] e)\s') = C F 5 ss , 



1 - z 



f*P _ 2- 



'1-Z k 2 



+ 



(C.13) 



(p-\Pgg{z,k ± ]P,n;e)\v) 



2C A 



1 - z 



au K ± K A 



+ z(l-z)g 



a/3 K j_ K ± 



l-z* k\ ' k\ 



+ 



+pr (d^kikl + g^k^kl - g av klk{ - g^k\k c A 



+ 



(C.14) 



Using the conventions for spinor products and polarization vectors of Ref. [41| and the rep- 
resentation of the transverse momenta in Eq. ( |6.56| ), we can derive the following relations, 
valid in the collinear limit, 



<A|/z) <T HA') = -5ay , (A|/i) (X'\a)g^ = -<y A( _ v) , 



k^ k v 

2 (Mv) HA') = -2 (A|a) ^ </?|A') = 6 XX , + *A(-V)^p 



Aj I ft/ I 



and 



-2 <A|//> (-A|a) 



A./ I a 



1. 



(C.15) 



(C.16) 



(C.17) 



Using these relations when inserting Eqs. ( |C.13| ) and ( |C.14| ) into Eq. ( C.12| ), we obtain the 
same results as in Eqs. ( |C,8| ) and (p.9|) . 
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